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PREFACE 


1.  The  Corps  of  Engineers,  tlirongli  its  Civil  Works  [)rograni,  has  sjjoiisored,  over  the  past 
23  years,  research  into  the  behavior  and  cliaracteristics  of  tidal  inlets.  The  Corps’  interest  in 
tidal  inlet  research  stems  from  its  responsibilities  for  navigation,  beach  erosion  prevention 
and  control,  and  flood  control.  Tasked  with  the  creation  and  maintenance  of  navigable  L'.S. 
waterways,  the  Corps  routinely  dredges  millions  of  cubic  yards  of  material  each  year  from 
tidal  iidets  that  connect  the  ocean  with  bays,  estuaries,  and  lagoons.  Design  and 
construction  of  navigation  improvements  to  existing  tidal  inlets  are  an  important  part  of  the 
work  of  many  Corps’  offices.  In  some  cases,  design  and  construction  of  new  inlets  are 
re(iuired.  Development  of  information  concerning  the  hydraulic  characteristics  of  inlets  is 
important  not  only  for  navigation  and  inlet  stability,  hut  also  because  inlets,  by  allowing  for 
the  ingress  of  storm  surges  and  egress  of  flood  waters,  play  an  importatit  role  in  the  flushing 
of  hays  and  lagoons. 

2.  A research  program.  General  Investigation  of  Tidal  Iidets  (On’l),  was  developed  to 
provide  ({uanlitative  data  for  use  in  design  of  inlets  and  inlet  improvements.  It  is  designed  to 
meet  the  following  objectives: 

To  determine  the  effects  of  wave  action,  tidal  flow,  and  related  phenomena  on  iidet 
stability  and  on  the  hydraulic,  geometric,  and  sedimentary'  characteristics  of  tidal 
inlets;  to  develop  the  knowledge  necessary  to  design  effective  navigation 
improvements,  new  inlets,  and  sand  transfer  systems  at  existing  tidal  iidets;  to  eviduate 
the  water  transfer  and  flushing  capability  of  tidal  inlets;  and  to  define  the  jirocesses 
controlling  inlet  stability. 

3.  The  GITI  is  divided  into  three  major  study  areas:  (a)  inlet  classification,  (b)  inlet 
hydraulics,  and  (c)  inlet  dynamics. 

а.  Inlet  Classification.  The  objectives  of  the  inlet  classification  study  are  to  classify 
inlets  according  to  their  geometry,  hydraulics,  and  stability,  and  to  determine  the 
relationships  that  exist  among  tlie  geometric  and  dynamic  characteristics  and  the 
environmental  factors  that  control  these  characteristics.  The  classification  study  keeps  the 
general  investigation  closely  related  to  real  inlets  and  produces  an  important  inlet  data  base 
useful  in  documenting  the  characteristics  of  inlets. 

б.  Inlet  Hydraulics.  The  objectives  of  the  inlet  hydraulics  study  are  to  define  the 
tide-generated  How  regime  and  water  level  fluctuations  in  the  vicinity  of  coastal  inlets  and 
to  develop  techni(|ue8  for  predicting  these  phenomena.  The  inlet  hydraulics  study  is  divided 
into  three  areas:  (1)  idealized  inlet  model  study,  (2)  evaluation  of  state-of-the-art  physical 
and  numerical  models,  and  (3)  prototype  inlet  hydraulics. 
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(1)  Tlic  Idfalizcd  Iidct  Model.  Tlic  objectives  of  this  model  study  are  to  determine 
the  effeet  of  iidet  configurations  and  structures  on  discharge,  head  loss  and  velocity 
distrihulion  for  a numher  of  realistic  iidcl  shapes  and  tide,  conditions.  An  initial  set  of  tests 
in  a trapezoidal  iidet  was  conducted  between  1967  and  1970.  However,  in  order  that 
suhseipieiit  inlet  models  are  more  representative  of  real  inlets,  a numher  of  “idi'alized” 
moilels  ri'presenting  various  inlet  morphological  classes  an;  being  develofied  and  hasted. 'I’lu' 
effects  of  jetties  and  wave  action  on  the  hydraulics  are  included  in  the  study. 

(2)  F.valuation  of  State-of-the-Art  Modeling  Tcchni(|ues.  The  objectives  of  this 
part  of  the  inlet  hydraulics  study  are  to  determine  the  usefulness  and  reliability  of  existing 
physical  and  numerical  moileling  techniipies  in  predicting  the  hydraulic  characteristics  of 
inlet-bay  systems,  and  to  determine  whether  simple  tests,  performed  rapidly  and 
economically,  are  useful  in  the  eviduation  of  projiosed  inlet  improvements.  Masonboro  inlet, 
N.  C.,  was  selected  as  the  prototype  inlet  which  would  be  used  along  with  hydraulic  and 
numerical  models  in  the  evaluation  of  existing  technicpios.  In  September  1969  a complete 
set  of  h\  (Iraulic  and  hafh\  metric  data  was  collected  at  .Masonboro  Inlet.  Construction  of  the 
fixed-bed  physical  model  was  initiated  in  1969,  and  extensive  tests  have  been  |)erformed 
since  then.  In  addition,  three  existing  numerical  models  were  applied  to  predict  the  inlet’s 
hydraulics.  Kxtensive  field  data  were  collected  at  Masonboro  Iidet  in  August  1974  for  use  in 
evaluating  the  capabilities  of  the  physical  and  numerical  models. 

(3)  Prototype  Inlet  Hydraulics.  Field  .studies  at  a number  of  inlets  are  providing 
information  on  prototx  pc  inlet-bay  tidal  hydraulic  relationships  and  the  effects  of  friction, 
waves,  tides,  and  inlet  inorjihologv  on  these  relationships. 

c.  Inlpl  Dynamics.  The  basic  objective  of  the  inlet  dynamics  study  is  to  investigate 
the  interactions  of  tidal  How,  inlet  configuration,  and  wave  action  at  tidal  inlets  as  a guide 
to  improvement  of  inlet  channels  and  nearby  shore  protection  works.  The  study  is 
subdivided  into  four  specific  areas:  (1)  model  materials  evaluation,  (2)  niovalilc-bed 
modeling  evaluation,  (3)  reanalysis  of  a previous  inlet  model  study,  and  (4)  prototype  inlet 
studies. 


(1)  Model  Materials  Evaluation.  This  evaluation  was  initiated  in  1969  to  provide 
data  on  the  resjionse  of  movable-bed  model  materials  to  waves  and  flow  to  allow  selection 
of  the  optimum  bed  materials  for  inlet  models. 

(2)  Movable-Bed  Model  Evaluation.  The  objective  of  this  study  is  to  evaluate  the 
state-of-the-art  of  modeling  teehniipies,  in  this  case  movable-bed  inlet  modeling.  Since,  in 
many  casi-s,  movable-bed  modeling  is  the  only  tool  available  for  [iredicting  the  response  of 
an  inlet  to  improvements,  the  capabilities  and  limitations  of  these  models  must  be 
established. 
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(3)  Reanalysis  of  an  Earlier  inlet  Model  Study.  In  1957,  a report  entitled, 
“Preliminary  Report:  Laboratory  Study  of  the  Effect  of  an  Uncontrolled  Inlet  on  tlie 
Adjacent  Beaches,”  was  published  by  the  Beach  Erosion  Board  (now  CERC).  A reanaly  sis  of 
the  original  data  is  being  performed  to  aid  in  planning  of  additional  GlTl  efforts. 

(4)  Prototype  Dynamics.  Field  and  office  studies  of  a number  of  inlets  are 
providing  information  on  tlie  effects  of  physical  forces  and  artificial  improvements  on  iidet 
morphology.  Of  particular  im|)ortance  are  studies  to  define  the  mechanisms  of  natural  sand 
bypassing  at  inlets,  the  response  of  inlet  navigation  channels  to  dredging  and  natural  forces, 
and  the  effects  of  inlets  on  adjacent  beaclies. 

4.  This  report  discusses  the  calibration,  base  tests,  and  predietive  tests  of  a numerical 
model  applied  to  Masonboro  Inlet,  N.  C.,  as  part  of  tbe  evaluation  of  tlie  state-of-the-art  of 
inlet  modeling  techniques.  It  presents  the  data  necessary  for  a comparison  of  results  of  the 
physical  and  numerical  models  discussed  in  the  basic  report  and  in  the  following 
appendixes: 

a.  Appendix  1.  Sager,  R.  A.,  and  Seabergh,  W.  C.,  “Fixed-Bed  Hydraulic  Model 
Results.” 

h.  Appendix  2.  Mascb,  F.  D.,  Brandes,  R.  J.,  and  Reagan,  J.  U.,  “Numerical 
Simulation  of  llydri'dynamics  (WRE)”  (In  2 Vols). 

c.  Appendix  3.  (3icn,  R.  M.,  and  Hembree,  L.  A.,  “Numerical  Simulation  of 
H\  drodynamics  (TR ACOR).” 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to  metric  (SI) 
units  as  follows: 


■Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

S(|uare  inches 

6.452 

stptare  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

s<]uare  feet 

0.0929 

square  meters 

eubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

s<|uare  yards 

0.836 

square  meters 

eubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

scpiare  miles 

259.0 

hectares 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

ciu  rices 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins* 

'To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings,  use  formula:  C = (5/9)  (F  — 32). 
To  obtain  Kelvin  (K)  readings,  use  formula;  K = (5/9)  (F  — 32)  +273.15. 
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SYMBOLS  AND  DEFINITIONS 


A 

^im 

"MSi, 


'.MTL 


\ 


il 

li. 


H. 


K, 

K,. 


variable  ba\  surface  area  (conslaiil  for  Keulegan  model) 

bay  surface  area  at  Iiigb  water 

bay  surface  area  at  low  water 

bay  surfaet'  area  at  mean  sea  level  datum 

bav  surface  area  at  mean  tide  level  datum 

averajre  surface  area  of  tidal  lagoon  or  bav 


a veloeitN  coefficient  defined  as  — 

2ir 

r 


XL\ . . 

H in  units  — 


semiranije  ocean  tide 

elevation  of  basin  water  surface  referred  to  MSL 
elevation  of  ocean  tide  referred  to  MSL 


s«-mirani;e  of  bav  tide 


time  derivative  of  bav  tide, 


dll, 

dt 


T a 


keulegan’s  dimensionless  coefficient  of  repletion,  . . , 

‘ 27rll  A-y/l  + LX/r 

eiiuivalent  Keulegan’s  coefficient  of  repletion  for  nonprismatie  inlets 


friction  constant  defined  as 


(1.4«6)^ 


L inlet  length 

M subscript  indicating  maximum  value 

Q total  inlet  discharge 

Qj  discharge  in  i*^*  channel 

T tidal  period  (usually  12.42  hours) 

V velocity  in  tidal  inlet  channel 
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SYMBOLS  AND  DEFINITIONS-Continued 

'a/, 

maximum  velocity  in  i'^'  channel 

velocity  as  function  of  x in  i*^  channel 

a 

tidal  inlet  cross-sectional  area 

inlet  area  at  minimum  cross  section 

‘‘mi 

minimum  cross-sectional  area  of  i^*  parallel  channel 

“r, 

area  as  function  of  x in  i*^  cliannel 

“o 

inlet  area  at  .MSL 

a 

o 

gravitational  acceleration 

11, 

dimensionless  bay  tide,  — 

'‘i  + 

floodtide,  maximum  bay  response,  dimensionless 

h,_ 

ebbtide,  maximum  bay  response,  dimensionless 

h, 

average  dimensioidess  bay  tide 

•>2 

"2 

dimensionless  ocean  tide,  — 

. 

channel  index  (i  = 1,  2,  3, . . . iy  ) 

hi 

the  maximum  number  of  cbannels 

j 

cross-sectional  index  (j  = 1 , 2,  3,  . . . ) 

maximum  number  of  cross  sections 

m 

subscript  indicating  minimum  value 

n 

Manning’s  roughness  coefficient 

‘1 

discharge  or  How  into  tidal  bay 

r 

inlet  hydraulic  radius 

effective  inlet  hydraulic  radius  at  MSL 

1 1 

SYMBOLS  AND  DEFINITIONS-Contiiiued 


li\  (Iraulic  radius  of  inlet  at  MSL 
r^j  liydraulie  radius  of  seetion  of  i</'  channel 

r^,  li\  (Iraulic  radius  as  function  of  x in  i'*  channel 

t time  from  iH'ginning;  of  tidal  cycle 

II  dimensioidess  factor  used  in  tiic  determination  of  inlet  velocity 

(u  =\/h2  — hj  ^ for  noodtide;u  = — x/ii,  “lij,  for  ebbtide) 

u+  maximum  value  of  u for  floodtide 

u—  maximum  negative  value  of  u for  ebbtide 

Wg  inlet  width  at  MSL 

X distance  measured  along  the  channel  axis 

27rt 

6 dimensionless  tidal  time, 

T 


TT  3.14159... 

2 sum  of  following  indexed  items 

n volume  of  tidal  prism  through  the  inlet,  2Hj^ 

a lag  of  bay  tide  behind  oeean  tide 

a+  tloodtide,  lag  of  bay  tide  bebind  oeean  tide 

a—  ebbtide,  lag  of  bay  tide  behind  ocean  tide 

|3  parameter  to  characterize  bay  surface  area  variation 

f inlet  width  variation  parameter,  inlet  beach  side  slope,  horizontal-vertical 

rj  average  tide  height  through  inlet 

X dimensionless  friction  factor  as  expressed  by  Keulegan, 

2gn^  _ 2gn^ 

(3.2808)V3  r'/3  “ (1.486)^  r‘/3 
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SIMPLIFIED  NUMERICAL  (LUMPED  PARAMETER)  SIMULATION 


by 

C. ,/.  Iliivul 
(unl 

(>.  L.  U iiilerfirrsl 
1.  INTRODUCTION 

This  study  is  cotiecrm'd  with  tlie  iiii|il('iiuMituti()ii  and  application  of  a liydraidic 
inatiu'inatical  tnodcl  for  predicting  ocean  tide-induced  current  velocities  within  a coastal 
iidet  and  the  wat('r  level  fluctuation  in  an  adjoining  embaytnent.  The  stud)  of  mathematical 
models  is  part  of  a more  general  inv«'sligation  of  tidal  inlets  conducted  by  the  U.S.  Army, 
(iorps  of  Engineers,  and  is  inl('nded  to  e\aluate  the  degree  to  which  mathematical  and 
physical  models  can  he  used  to  predict  (piatititativcly  the  hydrodynamics  (other  than 
sediment  moyement)  of  How  through  tiihd  inlets.  As  a typical  problem  area,  Masonhoro 
Inlet,  North  Carolina,  was  s<‘leclcd  as  a case  study  for  extensive  field,  mathematical  model, 
and  [)hysical  model  investigations. 

The  mathematical  mod<’l  used  in  this  study,  referred  to  as  the  lumjied  imrameter 
niiliroach,  is  based  on  an  extension  of  the  method  develo|jed  by  Keulegan  (1967).  Although 
the  method  reipiires  more  effort  to  im[)lement  and  apply  than  Keulegan ’s  method,  it  is 
relativeh  sim|ile  to  use  in  comparison  to  the  more  complete  schemes  whieh  are  presently 
available.  Moreover,  in  most  cases  it  should  give  a better  solution  of  inlet  velocities  and  basin 
responsi-  than  Keulegan's  method  sinee  it  ineor[)orates  more  of  the  physical  processes 
involved. 

The  numerical  system  described  in  this  study  is  composed  of  three  computer  programs, 
each  j)erforming  a separate  function.  One  program  generates  a set  of  tables  to  give 
generalized  iidet  hydraulics  for  some  variable  basin  surface  areas.  A second  [irogram 
(INLET)  gives  serial  calculations  of  the  inlet  How  and  the  basin  variations.  The  third 
|irogram  (SECPLT)  plots  the  ocean  tide,  basin  tidal  response,  inlet  velocity  and  inlet  flow 
and  computes  iidet  cross-sectional  areas  from  digitized  hydrographic  data. 

The  objective  of  this  study  is  to  apply  the  lumped  parameter  model  to  Masonhoro  Inlet 
and  determine  the  tidal  response  of  tlie  system  of  inner-connecting  channels  and  lelocities 
arising  from  a given  ocean  tide.  The  sco|)e  of  this  study  is: 

(a)  Mathematical  model  implementation  and  confirmation. 

(1)  Adapt  lumped  parameter  model  to  the  1969  bottom  survey  conditions. 

(2)  (ialibrate  the  model  for  prototype  currents  and  tides  of  12  September 

1969. 
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(It)  Matlicmatical  model  application. 

(1)  Predict  basin  tides,  tidal  [trism,  and  inlet  currents  for  jtreprojcct 
undeveloped  inlet  conditions  of  November  1964. 

(2)  Predict  basin  tides,  tidal  prism,  and  inlet  currents  for  modified  inlet  and 
north  jettv  conditions  of  July  1966. 

Ill  addition,  a description  of  the  model  development  is  presented  alon*;  with  .suggestions 
lor  a[)plying  the  model  to  other  inlets.  A documented  listing  of  the  computer  programs 
INLET  and  SECPLT  are  included. 

The  Masonboro  Inlet  study  area  (Fig.  1)  is  composed  of  an  entrance  channel  and  three 
inner-connecting  watercourses  referred  to  as  Masonboro  Channel,  Hanks  (Channel,  and  Shinn 
Oeek.  These  inner  channels  all  connect  with  the  Atlantic  Intracoastal  Waterwav.  The  mean 
high  water  (MIIW)  and  mean  low  water  (MLW)  lines  and  locations  where  velocitv 
measurements  were  made,  and  the  tidal  heights  observed  for  the  1969  survey  conditions  arc 
shown  in  Figure  1 . 

11.  BACKGROUND  ON  COMPUTATIONAL  METHOD 
The  subject  of  flows  and  corresponding  water  levels  of  tidal  inlet-bay  systems  has  been 
investigated  by  several  researchers.  Brown  (1928)  presented  a useful  solution  technique  for 
estimating  flows  tlirough  a coastal  inlet  and  the  variations  of  water  level  in  an  adjoining 
embay  ment.  Keulegan  (1967)  extended  this  work  and  gave  a more  useful  solution  based  on 
a dimensionless  form  of  the  equation  for  surface  changes.  Keulegan  offered  a simple  and 
eonvenient  tool  for  computing  the  main  hydraulics  of  tidal  inlets.  Because  the  method 
presented  in  this  study  is  based  on  an  extension  of  Keulegan’s  method,  a brief  review  of  his 
ba.sic  assumptions  and  pertinent  equations  is  given  in  this  section. 

An  inlet  system  with  a notation  consistent  with  Keulegan’s  formulation  is  shown  in 
Figure  2.  Kculegan’s  method  was  based  on  several  assumptions  related  primarily  to  the 
physical  characteristics  of  the  system  and  the  description  of  the  ocean  tide.  These  are 
identified  as: 

(a)  The  basin  responds  only  to  the  ocean  tide.  This  implies  that  the  basin  is 
bounded  except  at  the  inlet;  tlms,  there  is  no  freshwater  inflow  from  streams  or  lateral 
drainage  and  density  currents  are  negligible.  This  condition  also  implies  that  wind 
effects  are  neglected. 

(b)  The  ocean  tide  is  sinusoidal.  The  representation  of  ocean  tides  as  sinusoidal  is 
simple,  convenient,  and  often  reliable;  e.g.,  the  east  coast  tides  of  the  United  States  are 
well  represented  in  this  manner.  Along  the  U.S.  gulf  and  Pacific  coasts  and  in  many 
other  coastal  regions  of  the  world,  such  a representation  becomes  significantly  less 
accurate  because  the  tides  do  not  exhibit  simple  sinusoidal  form. 
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Kidiirc  1.  Masoiiljoro  Inlet  sliidv  area. 


(c)  Inertia  of  the  inlet  flow  is  negligible.  It  is  assumed  that  changes  of  flow  in  the 
inlet  ehannel  are  gradual  and  that  flow  aceeleration  can  be  neglected.  S|)ecifically,  the 
surface  slope  in  the  inlet  ehannel  is  balanced  by  the  velocity  head  and  channel 
resistanee. 

(d)  Length  and  width  of  bay  is  less  tlian  the  tidal  wavelength.  This  condition 
implies  that  the  basin  is  small,  and  if  deep  enough,  it  can  be  assumed  that  tlie  water 
surface  elevation  varies  uniformly  throughout  the  bay.  For  example,  if  a basin  is  scjuarc 
or  circular  (or  similar  in  form)  and  flow  is  not  st'riously  restricted,  there  may  he  a 
nearly  homogeneous  change  of  tlie  water  surface  elevation  over  the  entire  basin  in 
response  to  the  ocean  tide. 

(e)  The  inlet  is  prismatic.  The  cross  section  of  the  inlet  is  assumed  to  be  a uniform 
open  channel  with  entrance  and  outlet  well  defined.  However,  inlet  boundaries  are 
often  irregular. 

(f)  The  inlet  depth  is  large  eompared  to  the  tidal  range.  This  implies  that  the 
cross-sectional  area  of  the  inlet  is  considered  constant  over  the  tidal  cycle. 

(g)  Basin  walb  are  vertical.  This  implies  that  tlie  differential  storage  for  any 
differential  change  in  water  surfaee  elevation  within  the  basin  remains  eonstant. 

Based  on  the  above  assumptions,  it  is  possible  to  analyze  the  inlet-basin  response  for  a 
given  variation  in  the  ocean  tide.  Keulegan  (1967)  derived  the  following  relations  for 
estimating  the  surfaee  fluctuations  in  a basin  when  the  surfaee  of  the  sea  is  at  a higher 
elevation  or  a lower  elevation  than  the  surface  of  the  water  in  the  basin. 

For  floodtide: 

^ - KVh,  - 1,,',  h.  > I,,  . (1) 


For  ebbtide: 


±1 

Ad 


— K Vh|  — hj  , h|  > hj  , 


where,  K is  Keulegan’s  dimensionless  coefficient  of  repletion  and  defined  by: 

K = / 2gH 

27rH  A + LX^  ’ 


(2) 


(3) 
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II 


a 


I 

L 

A 


(lifiH'iisionle.s.s  ba\  tide  = , 

II 

(linu'tisionifss  ocean  tide  = , 

H 

27rl 

T ’ 

scniirani;c  ol'  ocean  tide  , 

tidal  period:  1 2.42  luiurs  for  tide  , 

2gn^ 

dimensionless  friction  coefficient:  r . , 

(1.486)2r'/^ 

vvliere  r is  measured  in  feet,  and  g is  measured  in  feet  per  square 
second.  Factor  1.486  is  tlie  conversion  from  metric  units  and  has 
the  dimensionsplengtli)'/3/tim(^(Cliow,  1959)  . 

inlet  hydraulic  radius  , 

Manning’s  resistance  coefficient  , 

inlet  cross-sectional  area  , 

time  from  heginning  of  tidal  cycle  , 

inlet  length,  and 

hay  or  lagoon  surface  area  . 


An  analv  tic  solution  to  the  above  eijuations  was  found  by  a Fourier  series.  Tables  by 
Keulegan  (1967),  give  ri-lations  to  predict  the  bay  ma.ximum  amplitude  (equal  for  floodtide 
and  ebbtide),  maximum  inlet  velocity  (equal  for  floodtide  and  ebbtide),  and  the  bay  tidal 
prism. 


m.  VARIABLE  BAY  SURFACE  AREA 

Van  de  Kreeke  (1967)  |)resented  a numerical  solution  of  the  inlet-bay  problem  which 
incorporated  the  effects  of  freshwater  inflow.  Mota  Oliveira  (1970)  formulated  and  solved 
numerically,  an  inlet  problem  which  accounted  for  variable  inlet  depth  and  variable  bay 
surface  area.  .A  numerical  computation  method  using  bay  inflow  with  solution  curves  was 
reported  by  Clenne,  Goodwin,  and  Glanzman  (1971).  Mota  Oliveira  (1970)  indicated  the 
importance  of  including  the  effect  of  varialde  bay  surface  area  in  the  computation  because  it 
can  have  a significant  effect  on  bay  tidal  responsi'  and  inlet  velocity.  He  determined  the  bay 
surface  area  as  a linear  function  of  the  bay  surface  elevation,  and  the  inlet  cross-sectional 
area  by  assuming  that  the  inlet  water  surface  elevation  is  the  average  of  ocean  and  bay 
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elevations.  Because  inlet  depth  was  allowed  to  vary  it  was  possible  to  vary  the  inlet  area  with 
respect  to  time.  It  was  shown  that  inclusion  of  these  improvements  in  the  computational 
scheme  would  give  better  results  for  predicting  inlet  hydraulics. 

Although  Keulegan  (1967)  did  not  account  for  a variation  of  bay  surface  area,  his 
method  is  useful  because  of  convenience.  The  collection  of  inlet  parameters  of  a single 
coefficient,  the  coefficient  of  repletion,  allows  solution  values  for  the  governing  equation  to 
be  conveniently  tabulated.  On  this  basis,  Keulegan’s  method  is  used  here  but  extended  by 
combining  Mota  Oliveira’s  scheme  for  including  the  variation  of  bay  surface  area.  Such  a 
solution  technique  would  retain  the  essential  simplicity  of  Keulegan’s  original  solution  but 
reduce  tlie  effect  of  the  most  restrictive  assumption,  that  of  vertical  bay  walls. 

Figure  3 illustrates  the  hnear  variation  of  bay  surface  area  ^v^th  surface  elevation  of  the 
bay  as  suggested  by  Mota  Oliveira  (1970).  The  varying  surface  area  is  given  by  ; 

A = A„  (1 +I3h,'  , (4) 

where  |3  is  the  hnear  bay  surface  area  variation  parameter,  and  A^  is  the  average  surface 
area  in  tlie  bay. 
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If  higli  and  low  water  areas  can  be  determined  from  hydrographic  cliarts  and  topograpliic 
maps,  may  be  taken  as  the  mean  surface  area.  Moreover,  the  surface  area  variation 
parameter,  p,  can  then  be  estimated.  Tlie  area-water  level  relationship  is  often  nonlinear  for 
many  problems  encountered  in  practice.  If  the  volume  of  the  tidal  prism,  fl,  has  been 
determined  for  the  system  from  discharge  measurements  at  the  inlet  and  the  tidal  range  in 
the  bay  is  also  known,  then  the  value  of  can  be  computed  from: 


A 


o 


n 


2H 


im 


(5) 


where,  is  the  tidal  semirange  in  tlie  bay.  A corresponding  value  for  P should  be 

obtained  from  the  best  available  information  relating  surface  area  A and  bay  elevation  Hj . 

To  include  the  effects  of  the  variation  of  the  bay  area,  a direct  substitution  of  eejuation 
(4)  is  made  into  equation  (1).  This  gives  for  floodtide: 


_ T / xlf. T 

dd  2nH  A„(l+^h,)  yji  + Xy^  ^ ^ 


(6) 


Rearranging  terms  yield  the  simpler  expression: 

dhj  Vh  — h ' 
dd  ” (1  -l-i3h,)  ’ 

where  the  coefficient  of  repletion  is  now  based  on  the  average  bay  surface  area  and  is  given 
by: 

K = I 

2,H  A„  • '<*> 

This  gives  a differential  equation  of  the  surface  changes  for  the  system  in  parameters  K and 
P.  The  friction  factor,  X,  expressed  in  terms  of  Manning’s  n,  and  the  hydraulic  radius  is 
given  by: 


2gn^ 

(1.486)^r'/3  • 


(9) 
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Eijuation  (7)  is  applicable  if  ^ h,  and  the  corresponding  equation  for  ebbtide 
or  h|  < hj  is: 


dh,  Vh,  -h,  ' 

de  (1  + -jh,)  ■ 


(10) 


Equations  (7)  and  (9)  have  been  solved  by  numerical  integration  for  various  values 
of  /3  and  k using  a simple  harmonic  ocean  tide  defined  by: 


bj  = sin  0 . 


(11) 


Tables  produced  from  these  solutions  are  presented  in  Appendix  A;  a typical  solution  with 
symbol  definition  is  shown  in  Figure  4.  The  use  of  the  tables  is  described  below. 

For  a well  defined  inlet,  the  determination  of  a,  r,  and  L is  relatively  simple;  n may 
be  estimated  from  tables  Manning’s  n and  from  the  available  information  about  bottom 
roughness  (more  general  cases  will  be  discussed  later).  These  data  can  be  used  with  tidal  data 
and  bay  area  characteristics  to  determine  K as  defined  in  equation  (8);  for  speeific  values 
of  K and  )3,  values  of  hj^.,  a+,  u+,  hj_,  a—,  u— , and  hj  can  be  found  from  the 
tables  in  Appendix  A.  The  first  six  values  provide  data  on  maximum  conditions  for  both 
floodtide  and  ebbtide.  For  /3  greater  than  zero,  h,  will,  in  general  be  nonzero. 
All  hj  values  can  be  dimensionalized  by  multiplying  the  dimensionless  form  by  the 
semirange  of  tide  in  the  sea,  H;  e.g.,  the  maximum  bay  response  is: 

= . (12) 

The  phase  lag  in  degrees  between  maximum  tide  in  the  sea  and  the  maximum  water  level  in 
tlie  basin  is  denoted  by  a (Fig.  4),  and  may  be  eonverted  to  time  by: 


At 


T 

a . 

360 


(13) 


Maximum  dimensionless  inlet  velocities,  u+  and  u— , at  floodtide  and  ebbtide,  are 
determined  from: 


(1  + 0h,)dh, 

ic  de" 


(14) 
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Typical  solution  and  definition  of  terms  for  variable  bay  solutions. 


From  coiitiiiuitv. 


27rll  (ili|  u 

do  ~ A ' ' 


(I'l) 


wlicrt-  V = tid;d  inlet  velocity-  Solviiii:  for  \ <;ives: 

27rH  A dll, 

~ T a dO  ■ 


(16) 


This  mav  also  be  written  in  terms  of  the  dimensionless  channel  velocity  as; 


27rH 

V = — - Ku  , 

I a 


(17) 


with  /5h,  = 0.0  or  in  the  eijnivalent  form: 


V = 


u . 


(18) 


IV.  THEORY  OF  GENERALIZED  LUMPED  PARAMETER 
In  addition  to  improving  the  computational  method  by  accounting  for  the  effect  ot 
variable  surface  area  in  the  basin,  it  is  also  possible  to  include  inertia  effects  in  the  How 
through  the  inl<-t  without  .significantly  increa.sing  the  computational  effort.  Although  inertia 
effects  are  generally  small  compared  to  the  principal  effects  which  govern  the  flow  through 
the  inlet,  their  inclusion  docs  provide  a more  complete  description  of  the  physical  processes 

involved.  Sliemdin  and  Forney  (1970)  developed  a model  which  included  the  el  feet  of  flow 
inertia  in  the  inlet.  The  lumped  parameter  model  developed  below  includes  both  the  inertia 
effects  and  .Mota  Oliveria’s  technique  for  including  the  effect  of  variable  surface  area  in  the 
basin. 

The  formulation  of  the  lumped  parameter  model  is  derived  here  from  energ>  and 
continuity  equations.  The  energy  equation  is  taken  as; 

Hj  - H,  = CJVi  V + ^ V , (19) 

where  heights,  II,  and  inlet  velocity,  V,  and  length,  L,  are  as  previously  defined, 
and  V represents  the  time  derivative  of  velocity  or  simply  dV/dt  (inlet  flow  acceleration). 
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A velocil)  cocl'ficiriit,  C^,,  introduced  for  notatioiial  simplicity,  accounts  for  inlet 
resistance  and  exit  losses.  This  coefficient  is  given  bv: 

C*  ^ 


C 


V 


(20) 


Based  on  procedures  suggested  by  Mota  Oliveira  (1970),  is  a function  of  the  head 
differential  (11^  — ll|)  since  the  hydraulic  radius,  r,  is  allowed  to  vary  with  time.  The 
value  of  r is  taken  approximately  as  the  sum  of  r^,  the  hydraulic  radius  at  mean  sea  level 
(MSL),  and  the  average  of  llj  and  ll^  . Introducing  r}  as  the  average  of  II j and  11^, 


V 


2 


(H,  +H,)  , 


(21) 


then 


= + 9 


(22) 


After  .substituting  for  X,  becomes: 


2g 


(1.486)2  ^ ° . 


(23) 


The  expression  of  conservation  of  mass  for  the  inlet-basin  system,  considering  that  the 
density  of  the  water  remains  constant,  is  given  by: 


AH|  = aV  + q , 


(24) 


where  the  overdot  signifies  the  time  derivative  and  (|  is  any  basin  inflow  which  does  not 
pass  through  the  inlet.  This  expression  is  an  equation  of  volume  continuity.  The  term  on  the 
left-hand  side  of  the  e(|uation  represents  the  change  in  volume  of  water  in  the  basin  during 
an  interval  of  time.  On  the  right-hand  side  of  the  eipiation,  the  first  term  is  the  discharge 
through  the  inlet  and  the  second  term  is  the  discharge  either  into  or  from  the  system  other 
than  through  the  inlet  during  the  same  time  interval.  Note  that  q is  positive  when  there  is 
inflow  to  the  system  and  negative  when  there  is  outflow  from  the  system.  Inflow  can  be 
caused  by  stream  discharge,  lateral  drainage  or  direct  rainfall  on  the  basin.  Outflow  may 
arise  from  water  being  discharged  from  the  basin  through  channels  which  lead  to  other 
bodies  of  water  of  lower  head.  Clearly  q is  to  be  taken  as  the  algebraic  sum  when  there  is 
both  inflow  and  outflow;  moreover,  ij  must  be  specified  as  a known  function  of  time. 
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E(juations  (19)  and  (24)  expressed  in  tlie  terms  of  tlie  time  derivatives  of  V and  Hj , are: 

V = ^ (H^ -H,  -C^  iVIV)  . (25) 

H.  = ^ (aV  + q)  . (26) 

Tliis  gives  two  simultaneous  differential  equations  in  which  both  and  (j  are  basic  input 
to  the  problem  and  prescribed  as  a function  of  time. 

Equation  (25)  is  nonlinear  because  energy  losses  in  the  inlet  channel  are  governed  by  the 
s(}uare  of  the  inlet  velocity  and  is  a function  of  the  head  differential  between  the  sea 
and  basin. 

-Mota  Oliveira  (1970)  assumed  vertical  walls  confining  tlie  inlet  which  yields  a 
cross-sectional  area  function  given  by: 

a = a^  -I-  w^T?  . (27) 

To  eliminate  tbe  restriction  of  vertical  inlet  walls,  a^  mav  be  taken  as  the  area  at 
MSL,  w^  is  taken  as  the  width  at  MSL  and  f is  introduced  as  the  beach  slope, 
horizontal-vertical,  within  the  tidal  range.  This  gives  a cross-sectional  area  function, 

^ • (28) 

It  has  been  found  that  for  beach  slopes  of  less  than  100:1,  the  quadratic  term  in  eijuation 
(28)  will  seldom  be  more  than  a few  percent.  Consequently,  it  is  concluded  that  a precise 
determination  of  f is  unnecessary. 

By  using  equations  (25)  and  (26)  both  the  bay  tide,  Hj , and  inlet  velocity,  V,  can  be 
estimated  if  the  following  functions  and  parameters  are  known: 


Hj(t), 

ocean  tide. 

q(t), 

inflow. 

mean  bay  surface  area. 

p, 

surface  area  variation  parameter. 

“o’ 

MSL  inlet  area, 

•■o’ 

MSL  inlet  hydraulic  radius. 

L, 

inlet  length. 

n, 

Manning’s  n. 

MSL  inlet  width,  and 

r. 

inlet  beach  slope. 
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As  lilt'll lioned  [iroviously , these  data,  otlier  than  Manning’s  n,  can  he  readily  determined 
fur  a well  defined  inlet;  their  deterriiination  for  the  more  common  irregularh  defined  inlets 
is  discussed  in  the  following  section.  Manning’s  n is  essentially  estimated  from  a knowledge 
of  hotlom  characteristics  and  later  adjusted  to  a satisfactory  calibration  of  the  model. 

V.  DETERMINATION  OF  EQUIVALENT  K FOR  NONPRISM ATIC  INLETS 

Recaiise  inlet  channels  are  generally  irregidar  in  both  widtii  and  dejith  throughout  their 
length,  the  assumption  of  a prismatic  inlet  is  not  always  valid.  To  account  for  such 
variations,  the  following  method  for  calculating  an  ecpiivalent  has  been  developed. 

The  inlet  vicinity  is  divided  into  several  parallel  channels  (three  shown  in  Fig.  5): 
partitioning  of  the  channel  is  arbitrarv.  The  goal  is  to  obtain  nearly  uniform  values  of  r,  X, 
V,  and  A in  each  cell,  but  generally  these  conditions  cannot  be  satisfied  exactly.  .\t  each 
cross  section.  A,  B,  C,  D . . . tile  cross-sectional  area  for  each  channel  1,  2,  3,  . . . is 
determined,  including  the  width  of  the  channel.  Along  each  shoreline  and  between  adjacent 
channels,  longitudinal  segment  length  can  be  measured  and  averaged  so  that  each  channel 
segment  is  given  a representative  length.  Since  an  area  and  width  are  known  for  each  section 
and  since  depths  arc  small  compared  to  widths,  representative  hydraulic  radii  can  be 
computed  for  each  section  from: 


area 

width 


(29) 


Table  1 illustrates  the  data  needed  in  defining  the  inlet  characteristics  for  the  inlet  shown 
schematically  in  Figure  5. 

The  coefficient  of  repletion  represents  the  effect  of  inlet  conditions  on  the  energy 
conversion  within  the  inlet.  This  eonversion  is  related  to  the  expression. 


(30) 


Within  the  parentheses,  the  first  term  represents  the  loss  of  energy  necessary  to  convert 
potential  energy  in  the  sea  to  kinetic  energy  in  tlm  inlet  channel.  The  velocity,  V,  in  the 
inlet  reaches  a maximum  value  where  the  contraction  is  the  greatest  or  at  the  location  where 
a minimum  cross-seciional  area  exists  in  the  inlet.  In  general,  this  minimum  inlet 
cross-sectional  area,  a^ , equals  tlie  sum  of  tlie  minimum  cross-sectional  areas  of 
the  i*^  channels, 

'm 

~ ^ i ’ 

i = 1 
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Figure  5.  Nonprismatic  inlet  representation. 


Table  1.  Required  data  specifications  for  nonprismatic  inlet  computation. 
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where  = luimber  of  ehannels.  For  the  example  in  Figure  5,  tlie  number  of  channels 
eijuals  three,  i^  = 3,  and  minimum  cross-sectional  areas  for  each  i*^  channel  are  at 
section  C,  so  a^-  = a^i.  For  this  example,  the  minimum  inlet  cross-sectional 
area,  a^,  is  found  by: 


a 


m 


3 

I = 1 


The  term  XL/r  represents  the  energy  lost  as  a result  of  bottom  friction.  Its  value  can  be 
determined  by  integrating  the  differential  effects  over  the  entire  inlet.  The  segmatized  inlet 
representation  facilitates  a numerical  approximation  to  such  a procedure. 

According  to  Keulegan  (1967),  the  energy'  loss  in  the  inlet  is  equal  to  the  head  difference 
between  sea  and  basin.  This  is  given  by: 

V2  = 2gH(h, -b,)  . (31) 

For  the  i*^  irregular  parallel  channel  this  may  be  rewritten  as: 


Mi 


0 


dx  = 2gH  (hj  - h, ) 


(32) 


where  V,  X,  and  r are  functions  of  x;  the  maximum  velocity  in  tlie  channel,  is 

related  to  the  flow  Q-.  Thus, 

Q,  = “mi  '>m  . 

and 


^*1  *ii  ^xi  • 

Substituting  equation  (33)  into  equation  (32)  gives: 


(34) 


Qi  = 


±+ 1 ^ 
aV  J r aV 


dx 


0 Yl  XI 


= 2gH  (h,  -h.) 


(35) 
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Solving  for  Q-  gives: 


= 


V2gH  (h^  - 1., )~ 


(36) 


Ki 


(lx 


r . a"”. 

0 Xt  XI 


By  combining  equations  (12)  and  (15)  the  continuity  ecjuation  may  be  expressed  in  the 
form: 


Q = 


(37) 


and  Q,  tlie  Bow  in  the  entire  inlet,  is  tlie  sum  of  the  flows  in  tlie  individual  channels, 


Also, 


Thus, 


dllj  _ 27rH  dhj 

dt  ~ T de 


27rHA  dh, 

T d^  ■ 


(38) 


(39) 


(40) 


Solving  e([uation  (40)  for  dhj/dg,  substituting  from  equation  (36)  and  multiplying  both 
numerator  and  denominator  by  a^  give: 


dll, 


a™  1 > 

f 

1 

~ 

i 

m 

a • 

2 Lx 

* I f' 

m 

a • 

2 

dx 

L ^ 

mi 

0 

XI 

K. 


(41) 


To  determine  K, , values  of  A,  a^,,  a^,,  r^,  are  taken  at  mean  tide  level.  The  integral  in 
e(|uation  (36)  is  replaced  by  the  numerical  approximation 


4 / ; 


XI  XI 


dx 


2^ 


1.486)^ 


■413 


2 

j = 1 


(Xy-Xy_,)  X„  = 0 


(42) 
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Thus,  an  effective  coefficient  of  re|)letion  can  he  (letcrniined  witli  the  foliovviiif:  e({uation; 


K 


1 


T 

27rII 


(43) 


If  an  inh't  is  not  [irismatic,  the  (ictemiination  of  L,  and  f is  more  difficult.  An 
e(|nivalent  coefficient  of  repletion  can  he  delermincd  as  described  in  tlic  [)revious  section. 
For  reasons  similar  to  those  discussed,  Wq,  and  f should  be  determined  in  the  throat 
repon:  a^^  being  the  MSL  iidct  area  at  the  throat,  \v^  the  MSI.  width,  and  f the 
representative  beach  slope  of  the  inlet  between  high  and  low  water. 

The  length,  L,  may  be  adetpiately  defined  as  the  average  length  of  the  parallel  channels 
selected.  Oidy  r^  remains  to  be  determined  and  the  following  procedure  may  be  used. 
Since  may  be  computed  from  ecjuation  (44),  r^^  may  be  taken  as  the  unknown  in  the 
original  definition  of  K (eq.  8). 

Substitution  from  (‘(piations  (8)  and  (9)  rearranging  terms  and  scpiaring  both  sides  of  the 
resulting  etpiation  gives: 


27rllAK. 


Ta., 


vW 


1 + 


(1.486) 


7 


4/3 


where  kj 


has  been  substituted  for  k in  equation  (8),  and  finally. 


1.104 

gn^L 

211  \ 

V 


- 1 


■3/4 


(44) 


Application  of  the  lumped  parameter  technicpie  is  discussed  in  Appendix  B. 

VI.  APPLICATION:  MASONBORO  INLET,  NORTH  CAROLINA 
The  calibration,  verification,  and  application  of  the  mathematical  model  described  to 
Ma.sonboro  lidet.  North  Carolina,  are  discussed  in  tliis  s»‘ction  (Fig.  6).  The  model  was 
calibrated  by  reproducing  tlie  prototype  inlet  velocities  of  12  September  1%9.  Apjdieation 
was  then  made  to  the  sy.stem  using  hydrographic  conditions  of  November  1964  and  July 
1966. 

Hydrographic  conditions  in  the  immediate  vieitiity  of  .Masonboro  lidet  in  October  1969 
are  sbown  in  Figure  7.  Tbe  grid  shown  in  this  figure  indicates  the  independent  channels 
adopted  for  representing  an  e(|uivalcnt  prismatic  inlet. 
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A tabulation  of  measured  prototype  velocities  for  12  September  1969  taken  at  the  inlet 
throat  is  included  in  Table  2.  These  velocities  were  used  with  the  throat  cross  section 
(Fig.  8)  to  compute  the  flows  given  in  Table  2 and  the  tidal  prism  in  Figure  9.  An  attempt 
was  made  to  compute  the  tidal  prism  from  the  bay  surface  areas  but  tl>is  method  did  not 
give  satisfactory  results.  This  is  probably  due  to  the  problem  of  delineating  the  effective 
surface  area  which  contributes  to  the  tidal  prism.  On  this  basis  tbe  inlet  velocity  data  were 
used  to  determine  the  tidal  prism.  A bay  tidal  range  of  4 feet  was  assumed  for  calculating 
the  mean  bay  surface  area  from  the  observed  ocean  tide  for  12  September  1969  shown  in 
Figure  10.  While  the  flows  shown  in  Figure  9 indicate  a small  outflow  through  other  outlets, 
the  magnitude  is  so  slight  in  comparison  to  the  total  flow  that  its  effect  is  considered 
negligible.  If  sufficient  information  was  available  for  describing  this  outflow  it  could  have 
been  included  in  the  numerical  computations. 

Using  the  procedure  described  previously,  tbe  grid  of  Figure  7 was  cross-sectioned  and 
segmented  so  that  the  equivalent  prismatic  inlet  characteristics  could  be  found.  The  cross 
sections  used,  indexing  from  ocean  to  lagoon,  are  shown  in  Figures  11  through  15.  These 
sections  were  |)lotted  and  their  areas  computed  by  the  program  SECPLT  (App.  C).  The 
various  areas,  widths,  and  lengths  determined  are  listed  in  Table  3 with  the  resulting 
e(|uivalent  area,  width,  length,  and  hydraulic  radius  produced  by  program  INLET  (Apj).  I)). 
The  values  thus  obtained  are  then  used  to  define  an  inlet  system  which  can  be  analyzed  bv 
the  subroutine  INLT2  whose  output  is  tabulated  in  Table  4.  An  initial  computation  was 
made  using  a simple  sinusoidal  tide  with  a period  of  12  hours  and  25  minutes.  A graphical 
presentation  of  this  analysis  is  shown  in  Figure  16  as  produced  by  subroutine  GRAPHIC; 
this  figure  permits  a comparison  of  computed  and  measured  prototype  velocities.  The 
relatively  good  velocity  agreement  was  obtained  by  [)ositioning  the  grid  in  Figure  7 so  that 
all  of  the  channels  at  any  given  cross  station  had  approximately  the  same  cross-sectional 
area  and  then  adjusting  the  value  assumed  for  Manning's  n to  obtain  maximum  agreement 
for  the  calibration  condition. 

In  a in‘C(ni(l  computation  for  the  12  September  1969  inlet  conditions  (Tables  5 and  6, 
and  Fig.  17),  the  obs<Tved  ocean  tide  was  used  as  input  rather  than  approximating  the  ocean 
tide  by  a .simple  sine  curve  of  appropriate  amplitude.  Overall  agreement  between  observed 
and  computed  velocities  was  not  substantially  changed.  Therefore,  in  later  computations  the 
ocean  tide  was  approximated  by  a sine  curve.  The  values  of  Manning's  n,  obtained  in  the 
calibration  procedure  were  retained. 

After  the  eomputations  for  12  September  1969  resulted  in  an  acceptable  comparison  of 
obsiTved  and  calculated  itdet  velocities,  the  lumped  ()arametcr  model  was  used  to  predict 
inlet-bay  interaction  for  various  conditions.  In  all  predictive  computations,  ecjuivalent  inlet 
dimensions  were  developed  by  the  same  procedure  ufk-d  with  the  12  September  1969  data 
and  the  same  Manning  c(M*fficient  (0.027)  was  used. 
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Table  2.  Prototype  data-velocities  and  flows  (Raii^e  2)  12  September  1%9. 


Measured  velocities  (ft/s) 


Time 

(h) 

Surface 

(hgt) 

Range  S 

Range  C 

Range  N 

Vel«M-it\ 

(meaii) 

Inflow 

M 

B 

M 

S 

m 

!!.,')() 

1.71! 

0.11 

1.0 

—0.5 

0.8 

1.0 

1.2 

1.8 

1.2 

1.2 

0 110 

15  111 

'>,(>() 

1 .41! 

1 

-1.0 

-1.6 

1.8 

-1.4 

-1.4 

-1.8 

1.2 

1.0 

1.0 

0.0.1 

15.12 

'(.fiO 

1.00 

- 2.5 

-3.5 

-3.5 

-2.8 

3.4 

-3.4 

0.8 

1.0 

1.0 

2.60 

11. (>0 

1 ().()() 

0.55 

-3.0 

-4.0 

-4.2 

-2.6 

-4.0 

-3.11 

-1.4 

-1.2 

-1.0 

-2.08 

45.75 

1 ()..■>() 

0.00 

-3.1 

-4.11 

- 4.0 

- 4.0 

-4,8 

-4.7 

0.8 

0.8 

-0.8 

-3.45 

-.50.4(> 

11.00 

-0.50 

-4.0 

—5.5 

—5.7 

-4.2 

-5,2 

—5.5 

0.4 

0.6 

0.8 

-3.80 

-53.20 

1 1 .,'>0 

-0.06 

-4.2 

-5.6 

-5.6 

-3.0 

-5.1. 

- 5.0 

-1.1 

-1.2 

0.6 

-3.83 

-51.25 

12.00 

-1.43 

-4.0 

-5.0 

-5.0 

-2.4 

-4.6 

-4.4 

-0.2 

-1.2 

-1.4 

-3.34 

42.62 

I2.r(0 

-1.112 

-4.3 

-4.0 

-4.9 

-2.8 

-4.2 

-4.2 

-1.0 

-0.4 

-0.2 

-3.18 

-38,98 

13.00 

-2.06 

-3.5 

-4.0 

-4.2 

-2.0 

—3.5 

-3.7 

-0.4 

-0.6 

- 0.6 

-2.61 

-31.22 

13.50 

-2.11 

-2.0 

-3.5 

-3.3 

—2.2 

-3.0 

-3.0 

-0.2 

-0.5 

-0.6 

-2.19 

-26.08 

14.00 

-2.10 

-2.4 

-2.6 

-2.4 

-1.8 

—2.2 

-2.2 

-0.5 

-0.3 

0.4 

-1.74 

-20.70 

14.50 

-1.94 

-1.0 

-1.4 

-1.5 

-1.2 

-1.2 

-1.2 

0.5 

0.4 

0.5 

-0.77 

-9.36 

15.00 

-1.54 

0.6 

0.6 

0.5 

0.5 

0.6 

0.5 

0.5 

0.6 

0.6 

0.57 

7.19 

I5..50 

-1.25 

0.6 

0.8 

1.2 

0.8 

1.2 

1.4 

1.0 

1.4 

1.4 

1 .08 

13.99 

10.00 

- 0.65 

1.3 

2.1 

2.0 

2.6 

2.6 

2.8 

1.8 

2.2 

2.8 

2.29 

31.61 

10.50 

0.00 

1.7 

2.6 

2.0 

2.8 

2.8 

2.6 

2.3 

2.7 

3.0 

2.57 

37.61 

17.00 

0..50 

2.5 

3.1 

3.4 

3.0 

3.6 

3.0 

2.5 

3.0 

3.3 

3.12 

47.77 

1 7. .50 

0.07 

2.4 

3.4 

2.6 

2.8 

3.6 

3.0 

2.6 

3.3 

3.6 

3.13 

49.86 

111.00 

l.iO 

2.6 

2.11 

2.6 

2.2 

3.8 

3.4 

2.2 

2.8 

3.3 

2.97 

48.53 

1»..50 

1.57 

2.11 

3.2 

3.0 

2.6 

3.4 

3.2 

1.8 

2.9 

3.4 

3.03 

,50.51 

10.00 

1 .00 

3.0 

.1.4 

3.0 

2.7 

3.4 

3.2 

2.0 

3.0 

3.2 

3.08 

52.74 

10.. 50 

2.10 

2.6 

2.11 

2.6 

2.0 

3.0 

3.0 

1.8 

2.5 

2.9 

2.66 

46.30 

20  00 

2 10 

2.0 

1.8 

2.2 

1.8 

2.0 

2.2 

1.8 

2.0 

2.0 

2.00 

34.69 

20.50 

1 .110 

0.6 

1.2 

1.3 

0.9 

0.8 

0.8 

0.7 

0.8 

1.0 

0.93 

15.91 

21. (Ml^ 

1.50 

-1.2 

-1.6 

1.11 

1 

1 

-1.2 

1 

2.0 

-1.8 

-0.7 

1 

-0,6 

1 

-0.8 

-1.39 

-23.09 

H ^ Ixjltorn 
VI  = rm*an  drptli 
S ==  siirfatT 
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Fiairr  H.  Katiw'  2 ijatiin};  t^lutioiis  and  assutncd  flow  s('cti()ns,  12  Srpli'tiilu  r l')6*F 


0 I 234  56789  10  II  12  0 

Time  ( h ) 


Figure  10.  Prototype  tidal  elevations,  12  September  1969. 
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Table  3.  Equivalent  coefficient  of  repletion  for  nonprismatic  inlet,  12  September  1%9. 


Verification  conditions 


Area  (ft^ ) 

Width  (ft) 

Length  (ft) 

9610 

9670 

6910 

8860 

1170 

1000 

800 

770 

0 

0 

0 

0 

5530 

3370 

5540 

6680 

940 

620 

760 

400 

1030 

1150 

1320 

1460 

4260 

6400 

5130 

4350 

820 

250 

265 

610 

860 

1180 

1400 

1570 

3860 

3900 

3850 

3030 

575 

165 

160 

410 

365 

585 

665 

725 

4020 

4020 

4020 

4020 

485 

160 

160 

340 

437 

532 

545 

550 

0 

0 

0 

0 

Tidal  period  = 

Tidal  semirange  = 

Bay  area  (mean  sea  level)  == 
Manning’s  n = 

Coefficient  of  repletion  = 
Inlet  area  = 

Inlet  length  = 

Inlet  depth  = 

Inlet  width  = 


12.417  hours. 

2.150  feet. 

1.928000E  + 08  square  feet. 
0.0270. 

1.553. 

1 .464000E  + 04  square  feet. 
3.593500E  + 03  feet. 

12.545  feet. 

1,310  feet. 


41 


Table  4.  Computed  tides  and  flows,  12  September  1969. 


Verification  conditions 


Inflow 

(M#f^s) 


HI 

(ft) 

Velocity 

(ft/s) 

Discharge 

(M#ft^8) 

-0.56 

3.10 

44.27 

-0.10 

3.38 

0.39 

3.52 

55.01 

0.55 

3.54* 

,56.01 

0.86 

3.50 

,56.98 

0.94 

3.48 

56.99* 

1.31 

3.28 

55.29 

1.70 

2.81 

48.71 

2.01 

2.02 

35.66 

2.05 

1.84 

32.68 

2.19 

0.71 

12.70 

2.21 ' 

0.05 

0.89 

2.15 

-1.54 

-27.14 

1.90 

-2.26 

-38.80 

1.61 

-2.64 

-43.88 

1.26 

-3.05 

-48.90 

0.85 

-3.39 

-52.07 

0.49 

-3.58 

-52.96* 

0.41 

-3.60 

-52.95 

0.01 

-3.67* 

-51.95 

-0.07 

-3.67 

-51.57 

-0.56 

-3.57 

-48.10 

-1.04 

-3.28 

-42.65 

-1.50 

-2.79 

-35.16 

-1.90 

-2.05 

-25.27 

-2.10 

-1.37 

-16.77 

-2.16 

-0.95 

-11.61 

-2.22* 

0.20 

2.45 

-2.18 

0.92 

11.21 

-1.95 

1.59 

19.80 

-1.66 

1.89 

24.24 

-1.32 

2.35 

31.12 

-0.93 

2.79 

38.42 

-0.49 

3.15 

45..38 

Critical  point  value. 

2.150  feel. 

1 .9280E  + 08  square  feet. 
0.350. 

0.03  feet. 

1.46400E  t 04  square  feet. 
1,310  feet. 

75:1. 

12.545  feel. 

3,594  feet. 

0.0270. 
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Tidal  sTmiraii^r  = 

Mean  bay  surface  area  = 

Bay  side  slope  parameter  beta  = 

Average  bay  level  = 

Inlet  propc'rties: 

X-»ection  area  below  mean  tide  level  = 
Inlet  width  at  mean  tide  level  = 

Inlet  beach  slope  = 

Mean  tide  level  depth  = 

Inlet  leiif^th  = 

Manning's  n = 


Table  5.  Kquivaleiil  coefficient  of  repletion  for  nonprismatic  inlet,  12  September  1969. 


Observed  tides 


Area(ft^) 

Width  (ft) 

Length  (ft) 

9610 

9670 

6910 

8tt60 

1170 

1000 

800 

770 

0 

0 

0 

0 

5530 

3370 

5.540 

6680 

940 

620 

760 

400 

1030 

11.50 

1320 

720 

4260 

6400 

5130 

4350 

820 

2.50 

265 

610 

860 

1 180 

1400 

75.> 

3860 

3900 

3850 

3030 

575 

165 

160 

410 

365 

.585 

665 

350 

4020 

4020 

4020 

4020 

485 

160 

160 

340 

437 

.532 

.545 

275 

0 

0 

0 

0 

Tidal  period  = 12.417  hours. 

Tidal  semirange  = 2.150  feet. 

Bav  ar<‘a  (mean  sea  level)  = I.92BOOK  + 04  square  feet. 
Manning's  n = 0.0270. 

(^oeffieient  of  repletion  = 1.636. 

Inlet  area  = 1.46400E  + 04  square  feet. 

Inlet  length  = 3.04250E  + 03  feel, 

lidel  depth  = 12.360  feet. 

Inlel  width  =1,310  feet. 


Table  6.  Computed  tides  and  Oows  for  nonprismatic  inlet,  12  September  1969. 


Observed 

tide  conditiom 

Time 

H2 

Inflow 

HI 

Velocity 

Discharge 

(h) 

(ft) 

(MiKt^s) 

(ft) 

(ft/«) 

(M#ft%) 

0.00 

0.13 

0.00 

-0.45 

3.37 

48.67 

0.63 

0.00 

0.04 

3.48 

52.55 

0.88 

0.00 

0.28 

3.52* 

54.32 

1.01 

0.00 

0.44 

3.49 

54.62* 

1.00 

1.07 

0.00 

0.52 

3.46 

54.35 

1.50 

1.36 

0.00 

0.96 

3.04 

49.38 

1.65 

1.94 

0.00 

0.00 

1.33 

1.67 

2.72 

2.58 

45.63 

44.47 

3.00 

2.08 

0.00 

1.96 

1.98 

34.79 

3.17 

2.09* 

0.00 

2.04 

1.61 

28.34 

3.50 

2.04 

0.00 

2.13 

0.56 

9.93 

3.67 

1.99 

0.00 

2.14* 

-0.24 

-4.23 

4.00 

1.87 

0.00 

2.05 

-1.74 

-30.41 

4.50 

1.61 

0.00 

1.81 

-2.07 

-35.47 

5.00 

1.27 

0.00 

1.54 

-2.36 

-.39.18 

5.50 

0.82 

0.00 

1.22 

-2.84 

-45.59 

6.00 

0.28 

0.00 

0.83 

-3.34 

-51.42 

-0.07 

-0.24 

0.00 

0.00 

0.54 

0.39 

-3.52 

-3.55 

-52.59* 

-52.32 

6.67 

-0.40 

0.00 

0.24 

-3.56* 

-51.68 

7.00 

-0.71 

0.00 

-0.07 

-3.53 

-49.92 

7.50 

-1.18 

0.00 

-0.55 

-3.44 

-46.73 

8.00 

-1.60 

0.00 

-1.03 

-3.26 

^2.52 

8.50 

-1.91 

0.00 

-1.49 

-2.81 

-35.52 

9.00 

-2.05 

0.00 

-1.88 

-l.% 

-24.22 

9.33 

-2.07* 

0.00 

-2.06 

-1.18 

-14.44 

9.50 

-2.07 

0.00 

-2.11 

-0.72 

-8.86 

9.75 

-2.05 

0.00 

-2.14* 

0.07 

0.81 

10.00 

-2.01 

0.00 

-2.10 

0.85 

10.37 

-1.68 

0.00 

-1.89 

1.70 

21.26 

10.83 

-1.40 

0.00 

-1.67 

2.18* 

27.96 

10.92 

-1.36 

0.00 

-1.61 

2.18 

28.05* 

11.00 

-1.32 

0.00 

-1.56 

2.14 

27.64 

-1.18 

-0.95 

0.00 

0.00 

-1.39 

-1.23 

1.98* 

2.13 

25.90* 

28.30 

12.00 

-0.34 

0.00 

-0.84 

3.00 

41.72 

12.50 

0.22 

o.tx' 

-0.37 

3.40 

49.43 

* Critical  point  value. 

Tidal  semirange 
Mean  bay  surface  area 
Bay  side  slope  parameter  beta 
Average  bay  level 
Inlet  properties: 

X-section  area  below  mean  tide  level 

Inlet  width  at  mean  tide  level 

Inlet  beach  slope 

Mean  tide  level  depth 

Inlet  length 

Manning's  n 


= 2.150  feet. 

= 1.928E  + 08  square  feet. 

= 0.350. 

= 0.05  feet. 

= 1.4640E  + 04  square  feet. 
= 1,310  feet, 

= 75:1. 

= 12.360  feel. 

= 3,042  feet. 

= 0.0270. 
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'I’lif  I'irfil  siK'li  ('onipiitution  was  for  loiig;-tcrin  mean  conditions  for  J 969.  Cross  sections 
for  1 2 SeplemlxT  1969  were  adjusted  by  0.35  foot  to  lower  the  reference  to  M.SL.  Then  llie 
mean  surface  area  of  tin'  bay  was  reduced  from  192,800,000  S(jnare  feet  computed  for  MTl. 
on  12  .''eplernber  1969  to  the  area  corresponding  to  MSL.  Tins  r<‘dnclion  was  obtained  by 
assumiim  that  the  area  varies  linearly  witli  the  water  level  and  is  t;overned  bv  a B of  0.35. 
'Hir  rneafi  surface  area  at  the  MSL  datum  is  given  by; 


A = A 

"MSL  "MTL 


2H 


= 1.866  X 10®  fl^  . 


where  b^  is  the  MTL  superelevation  of  0.35  foot  and  11  is  the  long-term  mean  semirange 
of  1.9  feel.  Other  iidet  characteristics  were  revised  accordingly  and  How  conditions 
computed  (Tables  7 and  8,  and  Fig.  18). 

Similar  eom|iutalions  were  made  for  the  following  conditions: 

(a)  1964  hydrography;  before  construction  of  north  jetty  (Figs.  19  and  20,  and 
Tables  9 and  10). 

(b)  1966  liydrograpliy;  after  north  jetty  was  constructed  (Figs.  21  and  22,  and 
Tables  1 1 and  1 2). 

(e)  Proposed  plan  H;  includes  construction  of  south  jetty  and  additioiuil  channel 
dredging  (Figs.  23  and  24,  and  Tables  13  and  14). 

With  all  conditions  modeled,  major  re.sults  (Table  15)  show  a general  reduction  in 
maximum  Hoodtide  velocities  from  the  1964  through  1969  test  conditions  (3.50  to  3.20 
feet  per  second).  The  introduction  of  the  proposed  south  jettv  would  bring  the  maximum 
Oijodtide  velocities  to  a higher  level  (3.33  IVel  per  second).  Maximum  ebbtide  velocities  are 
similarlv  affected. 

VII.  SUMMARY  AND  RECOMMENDATIONS 

A lumpeil  paranu'ter  mathematical  model  has  been  adapted  to  the  Masonboro  Inlet, 
North  Carolina,  system.  The  model  estimates  tidal  current  velocities  in  the  iidet  and  water 
level  fluctuations  in  the  adjoining  embay  menl.  f'alibrafion  and  verification  of  the  model  is 
carried  out  by  reproducing  prototype  conditions  of  12  September  1969.  Application  is 
made  for  |treprojecl  nndi’veloped  itdel  conditions  of  November  1964  and  modified  inlet  and 
north  jetty  conditions  id’ July  1966. 

The  lum[)ed  parameter  model  is  based  on  the  metliod  develojied  by  Keulegan  (1967),  but 
extended  to  include: 

(a)  V ariable  inlet  and  basin  surface  area, 

(b)  variable  inlet  channel  depths, 

(c)  variable  (nonsinusoidal)  ocean  tide, 

(d)  inlet  inertia  effects,  and 

(e)  bay  inllows  or  ontllows  other  than  the  sea. 
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Table  7.  Equivalent  coefficient  of  repletion  for  nonprismatic  inlet,  1969. 


Mean  sea  level  condition 


Area  (ft^) 

Width  (ft) 

Length  (ft) 

8350 

8290 

8460 

8280 

1080 

870 

1050 

725 

0 

0 

0 

0 

4990 

4240 

6710 

3950 

905 

995 

475 

330 

1030 

1150 

1320 

720 

3090 

6600 

4770 

4820 

750 

255 

250 

625 

860 

1180 

1400 

755 

4010 

4400 

3370 

2320 

550 

190 

150 

360 

365 

.585 

665 

3.50 

4550 

3220 

4510 

3280 

500 

130 

180 

345 

437 

532 

545 

275 

0 

0 

0 

0 

Tidal  period  = 

Tidal  semiraiige  = 

Bay  area  (mean  sea  level)  = 
Manning’s  n = 

Coefficient  of  repletion  = 
Inlet  area  = 

Inlet  length  == 

Inlei  depth  = 

Inlet  width  = 


12.417  hours. 

1.900  feet. 

1.86600E  + 08  square  feel. 
0.0270. 

1.681. 

I. 41000E  + 04  square  feet. 
3.04250E  + 03  feet. 

I I. -591  feet. 

1,250  feet. 
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Table  8.  Computed  tides  and  flows  for  nonprismatic  inlet,  1%9. 
Mean  sea  level  condition 


Time 

(h) 

H2 

(ft) 

Inflow 

HI 

(ft) 

Velocity 

(ft/s) 

Discharge 

0.00 

0.00 

-0.44 

2.82 

0.50 

0.48 

0.00 

3.07 

44.22 

1.00 

0.92 

0.00 

0.41 

3.20 

47.83 

1.17 

1.06 

0.00 

3.20' 

48.54 

1.50 

1.31 

0.00 

0.83 

3.16 

49.05' 

2.00 

1.61 

0.00 

1.22 

2.93 

46.97 

2.50 

1.81 

0.00 

1.56 

2.47 

40.50 

3.00 

1.90 

0.00 

1.82 

1.69 

28.18 

3.08 

1.90* 

0.00 

1.85 

1.52 

25.38 

3.50 

1.86 

0.00 

1.95 

0.38 

6.36 

3.58 

1.84 

0.00 

1.95* 

0.06 

0.96 

4.00 

1.71 

0.00 

1.87 

-1.58 

-26.16 

4.50 

1.45 

0.00 

1.65 

-2.03 

-32.87 

5.00 

1.09 

0.00 

1.39 

-2.41 

-37.93 

5.50 

0.67 

0.00 

1.07 

-2.80 

-42.75 

6.00 

0.20 

0.00 

0.71 

-3.11 

-45.70 

6.42 

-0.20 

0.00 

0.37 

-3.27 

-46.50* 

6.50 

-0.28 

0.00 

0.30 

-3.29 

-46.48 

6.83 

-0.59 

0.00 

0.01 

-3.33* 

-45.80 

7.00 

-0.74 

0.00 

-0.14 

-3.33 

-45.11 

7.50 

-1.16 

0.00 

-0.58 

-3.19 

-41.71 

8.00 

-1.50 

0.00 

-1.03 

-2.87 

-36.33 

8.50 

-1  74 

0.00 

-1.44 

-2.35 

-28.91 

-1.88 

0.00 

-1.77 

-1.58 

-19.03 

9.33 

-1.90' 

0.00 

-1.92 

-0.87 

-10.48 

9.50 

-1.89 

0.00 

-1.96 

-0.43 

-5.20 

9.67 

-1.87 

0.00 

-1.97' 

0.11 

1.31 

10.00 

-1.79 

0.00 

-1.89 

1.08 

12.% 

10.50 

-1.57 

0.00 

-1.68 

1.30 

15.92 

11.00 

-1.25 

0.00 

-1.44 

1.67 

20.% 

11.50 

-0.85 

0.00 

-1.13 

2.13 

27.55 

12.00 

-0.40 

0.00 

-0.77 

2.54 

34.01 

12.,50 

0.08 

0.00 

-0.37 

2.87 

.39.98 

^Critical  point  value. 

Tidal  semiraiige  =1.900  feel. 

Mean  bay  surface  area  = 1.866F^  + 08  s(]uarr  feet. 

Bay  side  slope  parameter  beta  = 0.350. 

Average  bay  level  = 0.02  feet. 

Inlet  properties: 

X-section  area  below  mean  tide  level  = 1 .4100K  + 04  wiuare  feet. 
Inlet  width  at  mean  tide  level  = 1,230  feet. 

Inlet  beach  slope  = 75:1. 

Mean  tide  level  depth  = 1 1. .591  feet. 

Inlet  length  = 3,042  feet. 

Manning's  n = 0.0270. 
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Kifjin'  19.  Hydrographic  survey  tor  November  1964. 
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condition,  1964. 


Table  9.  Etjuivalent  coefficient  of  repletion  for  nonprismatic  inlet,  1964. 

Mean  condition 


Area  (tV ) Width  (ft)  Length  (ft) 


5315 

10875 

8910 

10310 

1870 

1080 

1250 

1710 

0 

0 

0 

0 

4685 

4190 

4425 

3430 

1200 

420 

600 

960 

965 

1130 

1145 

545 

2890 

3095 

3330 

3250 

860 

185 

140 

340 

1075 

1375 

1375 

675 

5320 

4730 

4925 

5470 

1150 

310 

250 

520 

600 

800 

800 

400 

0 

0 

0 

0 

Tidal  period  = 12.417  hours. 

Tidal  STMiiirange  = 1.900  feet. 

Hay  area  (mean  sea  level)  = 1.86600E  + 08  square  feet. 
.Manning’s  n = 0.0270. 

Co»-fficient  of  repletion  = 1.520. 

Inlet  area  = 1.25650E  + 04  square  feet. 

Inlet  length  = 2.721 25E  + 03  feet. 

Inlet  depth  = 11.003  feet. 

Inlet  width  = 1,525  feet. 
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^Critical  point  value. 


Tidal  semirange 

Mean  bay  Hurfaee  area 

Hay  .side  slope  parameter  beta 

Average  bay  level 

Inlet  properties: 

X-seclion  area  b«*low  mean  tide  level 

Inlet  width  at  mean  tide  level 

Inlet  b<‘acli  slope 

Mean  tide  level  depth 

Inlet  length 

Manning's  n 


= 1.900  feet. 

= l.8()6K  + 08  s<iuare  feet. 

= 0.350. 

= 0.03  feel. 

= 1.2565E  + 04  si|uare  feet. 
= 1.525  feet. 

= 30:1. 

= 11.003  feet. 

= 2,721  feet. 

= 0.0270. 
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condition, 


Table  11.  Equivalent  coefficient  of  repletion  for  nonprismatic  inlet,  1%6. 


Mean  sea  level  condition 


Area  (ft* ) 

Width  (ft) 

* Length  (ft) 

7860 

6940 

6900 

9300 

1670 

940 

620 

665 

0 

0 

0 

0 

4790 

4840 

5170 

5060 

1210 

450 

370 

780 

1375 

1730 

1775 

880 

4450 

2630 

2970 

3510 

1050 

220 

150 

370 

1340 

1815 

1990 

1030 

5340 

6550 

6510 

4990 

1300 

400 

310 

440 

715 

1000 

1015 

505 

0 

0 

0 

0 

Tidal  period  = 

Tidal  semirange  = 

Bay  area  (mean  sea  level)  = 
Manning's  n = 

Coefficient  of  repletion  = 
Inlet  area  = 

Inlet  length  = 

Inlet  depth  = 

Inlet  width  = 


12.417  hours. 

1.900  feet. 

1.86600E  + 08  square  feet. 

0. 0270. 

1. -'iSS. 

1.35600E  + 04  stjuare  feet. 
3.79250E  + 03  feet. 

13.125  feet. 

1,790  feet. 
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Table  12.  Computed  tide  and  Oow  lor  nonpriamatic  inlet,  1966. 
Mean  sea  level  condition 


Time 

(h) 

H2 

(ft) 

Inflow 

(Mi(fft%) 

HI 

(ft) 

Velocity  j 

(ft/») 

Discharge 

(Mdft^s) 

0.00 

0.00 

-0.50 

2.93 

38.46 

0.48 

0.00 

3.17 

44.06 

1.00 

0.92 

0.00 

3.27 

48.19 

1.08 

0.99 

0.00 

3.27* 

48.67 

1.50 

1.31 

0.00 

0.78 

3.22 

49.93* 

2.00 

1.61 

0.00 

1.19 

2.98 

48.26 

2.50 

1.81 

0.00 

1.54 

2.51 

42.04 

3.00 

1.90 

0.00 

1.81 

1.74 

29.86 

3.08 

1.90* 

0.00 

1.85 

1.58 

27.10 

3.50 

1.86 

0.00 

1.96 

0.50 

8.63 

3.67 

1.82 

0.00 

1.97* 

-0.13 

-2.24 

4.00 

1.71 

0.00 

1.90 

-1.48 

-25.20 

4.50 

1.45 

0.00 

1.67 

-2.07 

-34.27 

5.00 

1.09 

0.00 

1.40 

-2.40 

-38.19 

5.50 

0.67 

0.00 

1.09 

-2.78 

-42.28 

6.00 

0.20 

0.00 

0.73 

-3.11 

-44.75 

6.33 

-0.12 

0.00 

0.47 

-3.26 

-45.25* 

6.50 

-0.28 

0.00 

0.34 

-3.32 

-45.16 

7.00 

-0.74 

0.00 

-3.40* 

^3.62 

7.50 

-1.16 

0.00 

-0.51 

-3.33 

-40.35 

8.00 

-1.50 

0.00 

-0.93 

-3.09 

-35.56 

8.50 

-1.74 

0.00 

-1.33 

-2.67 

-29.29 

9.00 

-1.88 

0.00 

-1.67 

-2.00 

-21.29 

9.33 

-1.90* 

0.00 

-1.84 

-1.39 

-14.57 

9.50 

-1.89 

0.00 

-1.91 

-1.01 

-10.53 

9.83 

-1.83 

0.00 

-1.96* 

0.02 

0.17 

10.00 

-1.79 

0.00 

-1.95 

0.72 

7.53 

10.50 

-1.57 

0.00 

-1.74 

1.61 

17.42 

11.00 

-1.25 

0.00 

-1.48 

1.87 

21.05 

11.50 

-0.85 

0.00 

-1.18 

2.27 

26.79 

12.00 

-0.40 

0.00 

-0.83 

2.66 

33.21 

12.50 

0.08 

0.00 

-0.43 

2.98 

39.46 

'critical  point  value. 

Tidal  semirange  = 1.900  feet. 

Mean  bay  surface  area  = 1.866E  + 08  square  feet. 

Bay  side  slope  parameter  beta  = 0.350. 

Average  bay  level  = 0.02  feet. 

Inlet  properties: 

X-section  area  below  mean  tide  level  = 1.3560E  + 04  square  feet. 
Inlet  width  at  mean  tide  level  = 1,790  feet. 

Inlet  beach  slope  =75:1. 

Mean  tide  level  depth  = 13.125  feet. 

Inlet  lengtii  = 3,792  feet. 

Manning's  n = 0.0270. 


Table  13.  Equivalent  coefficient  of  repletion  for  nonprismatic  inlet,  1969. 


PlanB 


Area  (ft^ ) 

Width  (ft) 

Length  (ft) 

3900 

3900 

3900 

3900 

485 

275 

235 

180 

0 

0 

0 

0 

5129 

5120 

5120 

5120 

585 

320 

465 

325 

1540 

1715 

1775 

900 

3430 

3430 

3430 

3430 

575 

145 

145 

415 

1860 

1890 

1925 

975 

4120 

4120 

4120 

4120 

570 

255 

170 

385 

615 

685 

680 

330 

0 

0 

0 

0 

Tidal  period  = 

Tidal  semirange  = 

Bay  area  (mean  sea  level)  = 
Mannings’  n = 

Coefficient  of  repletion  = 
Inlet  area  = 

Inlet  length  = 

Inlet  deptli  = 

Inlet  width  = 


12.417  hours. 

1.900  feet. 

1.86600E  + 08  square  feet. 
0.0270. 

1.802. 

1.37200E  + 04  square  feet. 
3.72250E  + 03  feet. 

16.703  feet. 

1,280  feet. 
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Table  14.  Computed  tide  and  flow  for  nonpriamatk  inlet,  1969. 


^Critical  point  value. 

Tidal  semirange  * 1.900  feel. 

Mean  bay  surface  area  = 1.866E  + 08  square  feet. 

Bay  side  slope  parameter  beta  = 0.350. 

Average  bay  level  = 0.02  feet. 

Inlet  properties: 

X-section  area  below  mean  tide  level  = 1.3720E  + 04  square  feet. 
Inlet  width  at  mean  tide  level  = 1,280  feet. 

Inlet  beach  slope  - 30:1. 

Mean  tide  level  depth  = 16.703  feet. 

Inlet  length  - 3,722  feet. 

Manning  s n 0.0270. 
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Table  15.  Summar)  of  conditions. 


Modeled  and  major  results 


Data 

Tides,  12  Sept.  1969 

Lung  term 

N.S.  jetties 

Sinusoidal 

Observed 

1969 

1966 

1964 

Plan  B 

II  ft 

•1(1) 

2.15 

2.15 

1.90 

1.90 

1.90 

1 .90 

192.8 

192.8 

186.6 

186.6 

186.6 

186.6 

A„Mft^ 

A' 

0.35 

0.35 

0.35 

0.35 

0.35 

0.35 

a„  ft  2 

14,640.00 

14,640.00 

14,100.00 

13,.560.00 

12,570.00 

13,720.00 

••o 

12.4 

12.4 

11.6 

13.1 

11.0 

16.7 

Lft 

3,040.00 

3,040.00 

3,040.00 

3,790.00 

3,720.00 

3,720.00 

W„  ft 

1,310.00 

1,310.00 

1,2.50.00 

1,790.00 

1,525.00 

1 ,280.00 

r 

75:1 

75:1 

75:1 

75:1 

30:1 

30:1 

Response ' 

II,  + ft 

2.14 

2.21 

1.95 

1.97 

1.94 

1.98 

II,-  ft 

2.14 

2.22 

1.97 

1.96 

1.93 

1.99 

V + ft/s 

3.52 

3.58 

3.20 

3.27 

3..50 

3.33 

V - ft/s 

3.56 

3.71 

3.33 

3.40 

3.61 

3..52 

Q + M#ft^/s 

54.6 

57.6 

49.1 

49.9 

49.0 

49.6 

Q - M#ft'Vs 

52.6 

57.6 

46.5 

45.2 

44.6 

47.8 

a + h 

0.5 

0.5 

0.5 

0.6 

0.6 

0.5 

a _ li 

0.4 

0.3 

0.3 

0.5 

0.5 

0.35 

'Maximum  values  + floudtide  — ebbtide. 
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The  method  is  based  on  a tidal  system  eomposed  of  a bay  eonneeted  to  tlie  open  oeean  by 
an  inlet  channel.  It  is  presumed  that  the  current  velocities  can  only  exist  in  the  inlet  channel 
and  tliose  in  the  ocean  and  bay  are  small  and  may  be  neglected.  Thus,  a rise  of  water 
elevation  in  the  sea  due  to  a tidal  action  results  in  a conversion  of  potential  to  kinetic  energy 
in  tlie  inlet  channel  and  tlien  part  of  the  kinetic  energy  is  converted  to  potential  energy  in 
the  bay  and  part  is  lost  due  to  energy  dissipation  in  the  channel.  In  the  bay,  this  means  that 
the  water  surface  varies  uniformly  over  the  entire  basin. 

The  lumped  parameter  mathematical  model  has  been  programed  for  a GE  400  computer. 
Two  digital  computer  programs,  one  for  tlie  numerical  computations  of  tlie  inlet-basin 
problem  and  tlie  other  for  plotting  and  computing  inlet  cross-sectional  areas,  are  provided 
together  with  documentation.  These  programs  sliould  be  generally  adaptable  to  otlier 
inlet-basin  problems  with  a minimum  of  effort. 

The  lumped  parameter  model  should  be  applicable  for  a variety  of  inlet  problems  where 
ocean  tides  penetrate  a tidal  bay  or  lagoon  through  a single  constricted  inlet.  The  model  is 
sufficiently  general  to  accommodate  a variety  of  inlet  shapes.  However,  care  should  be 
exercised  in  using  the  model  to  predict  water  motions  for  all  inlet-bay  problems.  Water 
motions  must  be  restricted  to  those  induced  only  by  ocean  tides  and  the  physical 
characteristics  of  the  system  must  approximately  conform  to  those  inferred  by  the 
underlying  assumptions  used  in  developing  the  model. 

It  is  believed  that  the  lumped  parameter  model,  while  generally  useful  for  estimating 
current  velocities  in  inlet  channels  and  water  level  fluctuations  in  a connecting  basin,  could 
be  improved.  The  introduction  of  two  coefficients  of  repletion,  one  for  floodtide  and  one 
for  ebbtide,  and  a different  computational  method  for  the  nonprismatic  inlets  to  inclu'de  the 
large  differences  in  floodtide  and  ebbtide  flows  tlirough  the  inlet,  are  possible 
improvements. 
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APPENDIX  A 


TABLES  OF  NUMERICAL  SOLUTIONS  FOR  VARIABLE  BAY  SURFACE  AREAS 


Table  1.  Numerical  solution  for  variable  bay  surface  areas  (Beta  = 0. 10). 


Repletion 

coefficient 

Flood  tide 

Ebbtide 

Average 

K 

H,' 

a 

U" 

H, 

a 

u 

Hi 

76.0 

0.987 

-0.237 

■■ 

-0.988 

0.001 

0.30 

■■ 

-0.974 

0.001 

0.40 

63.0 

mm 

-0.954 

0.002 

0.50 

0.546 

58.0 

WSSM 

■■ 

-0.929 

0.003 

0.60 

0.631 

50.0 

-0.646 

50.0 

-0.900 

0.003 

0.70 

0.706 

45.0 

0.852 

-0.724 

43.0 

-0.868 

0.004 

0.769 

40.0 

-0.790 

38.0 

-0.833 

0.005 

34.0 

-0.844 

32.0 

-0.797 

0.005 

1.00 

31.0 

28.0 

-0.760 

0.005 

1.10 

0.898 

26.0 

23.0 

-0.724 

0.005 

1.20 

0.924 

23.0 

0.668 

19.0 

-0.689 

0.005 

1.30 

0.945 

19.0 

0.635 

-0.964 

15.0 

-0.656 

0.004 

1.40 

0.960 

16.0 

0.604 

-0.977 

13.0 

-0.624 

0.004 

1.50 

0.972 

14.0 

0.576 

-0.986 

10.0 

-0.595 

0.003 

1.60 

0.981 

12.0 

0.550 

-0.992 

7.0 

-0.567 

0.003 

1.70 

0.987 

9.0 

0.526 

-0.996 

5.0 

-0.541 

0.002 

1.80 

0.991 

7.0 

0.503 

-0.998 

4.0 

-0.518 

0.002 

1.90 

0.995 

6.0 

0.482 

-0.999 

3.0 

-0.495 

0.002 

0.997 

5.0 

0.463 

-1.000 

1.0 

-0.475 

0.001 

2.10 

0.998 

4.0 

0.444 

-1.000 

1.0 

-0.456 

0.001 

2.20 

0.999 

2.0 

0.427 

-1.000 

0.0 

-0.438 

0.001 

2.30 

1.000 

2.0 

0.411 

-1.000 

0.0 

-0.421 

0.001 

2.40 

1.000 

1.0 

0.397 

-1.000 

0.0 

-0.406 

0.001 

2.50 

1.000 

1.0 

0.383 

-1.000 

0.0 

-0.391 

0.001 

2.60 

1.000 

1.0 

0.369 

-1.000 

0.0 

-0.377 

0.001 

2.70 

1.000 

0.0 

0.357 

-1.000 

0.0 

-0.364 

0.000 

2.80 

1.000 

0.0 

0.345 

-1.000 

0.0 

-0.352 

0.000 

2.90 

1.000 

0.0 

0.334 

-1.000 

0.0 

-0.341 

0.000 

1.000 

0.0 

0.324 

-l.O(K) 

0.0 

-0.330 

0.000 

* Ratio  bay  amplitude  to  the  sea  amplitude. 
^Dimensionless  maximum  inlet  velocity. 
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Table  2.  Numerical  solution  for  variable  bay  surface  areas  (Beta  = 0.25). 


Repletion 

coefficient 

Flood  tide 

Ebbtide 

Average 

K 

H,' 

a 

u" 

H| 

a 

u 

H, 

-0.239 

76.0 

-0.990 

0.001 

-0.355 

68.0 

-0.977 

0.003 

0.40 

0.446 

0.943 

-0.465 

63.0 

-0.959 

0.004 

0.50 

0.913 

-0.567 

56.0 

-0.936 

0.007 

0.60 

52.0 

0.878 

-0.659 

49.0 

-0.910 

0.009 

0.70 

0.694 

47.0 

0.841 

-0.739 

43.0 

-0.880 

0.010 

0.756 

41.0 

0.802 

36.0 

-0.847 

0.012 

0.868 

36.0 

0.763 

-0.862 

31.0 

-0.813 

0.013 

1.00 

0.849 

32.0 

0.725 

25.0 

-0.777 

0.013 

1.10 

0.883 

28.0 

0.689 

-0.938 

21.0 

-0.742 

0.012 

0.969 

24.0 

0.656 

-0.961 

16.0 

-0.708 

0.012 

0.931 

0.624 

-0.977 

13.0 

-0.675 

0.011 

1.40 

0.948 

0.595 

-0.988 

9.0 

-0.643 

0.010 

1.50 

0.961 

16.0 

0.568 

-0.994 

6.0 

-0.613 

0.008 

1.60 

0.971 

14.0 

0.544 

-0.997 

4.0 

0.007 

1.70 

0.979 

12.0 

-0.999 

3.0 

0.006 

1.80 

0.985 

10.0 

-1.000 

1.0 

0.005 

1.90 

0.989 

8.0 

0.479 

-1.000 

0.0 

0.004 

2.00 

0.992 

7.0 

0.461 

-1.000 

0.0 

0.004 

2.10 

0.995 

6.0 

0.443 

-1.000 

0.0 

-0.471 

0.003 

2.20 

0.997 

5.0 

0.427 

-1.000 

0.0 

-0.453 

0.003 

2.30 

0.998 

4.0 

0.411 

-1.000 

0.0 

-0.435 

0.002 

2.40 

0.999 

3.0 

0.397 

-1.000 

0.0 

-0.419 

0.002 

2.50 

0.999 

2.0 

0.383 

-1.000 

0.0 

-0.404 

0.002 

2.60 

1.000 

2.0 

0.371 

-1.000 

0.0 

-0.390 

0.001 

2.70 

1.000 

1.0 

-1.000 

0.0 

-0.376 

0.001 

2.80 

1.000 

1.0 

-1.000 

0.0 

0.001 

2.90 

1.000 

1.0 

-1.000 

0.0 

0.001 

3.00 

1.000 

0.0 

-1.000 

0.0 

-0.341 

0.001 

’Ratio  bay  amplitude  to  the  sea  amplitude. 
^Dimensionless  maximum  inlet  velocity. 
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Table  3.  Numerical  solution  for  variable  bay  surface  areas  (Beta  = 0.50). 


Repletion 

coefficient 

Flood  tide 

Ebbtide 

Average 

K 

Ht' 

a 

Hi 

a 

u 

H, 

0.232 

77.0 

0.983 

-0.242 

0.003 

0.339 

-0.362 

0.006 

0.40 

0.439 

65.0 

0.935 

-0.478 

-0.967 

0.009 

0.529 

58.0 

0.901 

-0.588 

54.0 

-0.948 

0.014 

0.609 

52.0 

-0.687 

47.0 

-0.925 

0.018 

0.70 

0.679 

47.0 

-0.773 

40.0 

-0.899 

0.021 

0.737 

43.0 

0.783 

-0.845 

32.0 

0.024 

0.786 

38.0 

0.744 

25.0 

-0.838 

0.025 

1.00 

0.826 

34.0 

0.709 

-0.941 

20.0 

-0.806 

0.025 

1.10 

0.859 

31.0 

0.676 

-0.969 

14.0 

-0.773 

0.024 

0.856 

28.0 

0.645 

-0.986 

10.0 

0.022 

25.0 

0.617 

-0.996 

5.0 

0.020 

1.40 

0.926 

23.0 

0.590 

-0.999 

2.0 

-0.678 

0.018 

1.50 

0.941 

20.0 

0.566 

0.0 

-0.649 

0.016 

1.60 

0.953 

18.0 

0.543 

-1.000 

0.0 

-0.621 

0.013 

1.70 

0.963 

15.0 

0.522 

-1.000 

0.0 

-0.594 

0.012 

1.80 

0.971 

14.0 

0.502 

-1.000 

0.0 

-0.569 

0.010 

1.90 

0.977 

13.0 

0.484 

-1.000 

0.0 

-0.546 

0.009 

0.982 

11.0 

0.467 

-1.000 

0.0 

-0.524 

0.007 

2.10 

0.986 

10.0 

0.450 

-1.000 

0.0 

-0.504 

0.006 

2.20 

8.0 

0.435 

-1.000 

0.0 

-0.485 

0.005 

2.30 

0.992 

7.0 

0.420 

-1.000 

0.0 

-0.467 

0.005 

2.40 

0.994 

6.0 

0.407 

-1.000 

0.0 

-0.450 

0.004 

2.50 

0.996 

5.0 

0.394 

-1.000 

0.0 

-0.434 

0.003 

2.60 

0.997 

4.0 

-1.000 

0.0 

-0.419 

0.003 

2.70 

0.998 

4.0 

-1.000 

0.0 

-0.404 

0.003 

2.80 

0.999 

3.0 

-1.000 

1.0 

-0.391 

0.002 

2.90 

0.999 

2.0 

-1.000 

1.0 

-0.379 

0.002 

3.00 

0.999 

2.0 

0.339 

-1.000 

0.0 

-0.367 

0.002 

* Ratio  bay  amplitude  to  the  sea  amplitude. 
2 

Dimensionless  maximum  inlet  velocity. 
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APPENDIX  B 


APPLICATION  OF  THE  LUMPED  PARAMETER  TECHNIQUE 

The  fuiiclamental  equations  of  the  lumped  parameter  approach  to  the  study  of  tidal  inlet 
hydraulics  arc  summarized  below: 


dV 

dt 

= -cjviv)  , 

(25b) 

dH 

1 , 

dt 

= ^(aV  + q)  , 

(26b) 

A 

= A„(l+^b,)  , 

(4b) 

c. 

= ^ll  + FL[r<,  + |(H,+H,  )]-*/.)  , 

(23b) 

F = 


2gn^ 


r„  = 


K.  = 


(1.486)2 

gn2 

1.104 

1 

g 

I--. 

2H 

T 

Ta„ 


- 1 


-3/4 


27rH  A 


i = 1 


Mm  ^ 

< 

L 

r..-V3  1 

Mmi  ) 

j=  1 

•V2 


(44b) 


(43b) 
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Numerical  constants  liavc  Im'cii  simplified  and  auxiliar)  values  introduced  to  simplifv 
programing,  but  each  eejuation  is  algebraically  e<|uivalent  to  an  e<|uation  already  identified 
in  the  text  by  the  same  arabic  numeral. 

Note  that  F is  changed  only  by  changes  in  .Manning’s  n;  Ki  depetuls  only  on  clianges 
in  F and  changes  in  the  manner  of  partitioning  the  itd«‘t  to  account  for  irrcgidar  cross 
sections.  Since  the  friction  term,  which  involves  F",  is  generally  smaller  than  tlic  otlier  terms 
in  S(|uare  brackets  in  the  definition  of  Kj  ; K, , itself,  is  not  very  sensitive  to  changes  iti  the 
manner  of  partitioning  the  inlet  or  to  changes  in  n.  The  effective  hydraulic  radius,  r,^.  is  a 
function  of  both  F"  and  Kj.  Since  F'  modifies  all  terms  in  r^^,  not  just  the  snudler  term, 
as  for  Ki,  r^  is  more  sensitive  than  Ki  to  changes  in  Manning’s  n.  F',  K],  and  r^^  are  all 
constant  for  any  one  value  of  Manning’s  n and  one  partitioning  sy  stem. 

The  calculation  of  re(|uires  F and  r^,  both  of  which  depcfid  on  Manning’s  n,  and 
the  tide  levels  in  both  bay  and  ocean,  ll|  and  II2.  The  tide  levels,  II,  and  ll^  ar<‘ both 
functions  of  time.  Thus  must  be  constantly  u|)dated  throughout  the  calculation. 

The  differential  c(juation  (25b)  is  nonlinear  because  is  a function  of 
both  11,  and  II2  and  because  a (juadratic  term  in  V a})pears.  The  differential  e(iuation 
(26b)  is  nonlinear  because  the  bay  area.  A,  is  a function  of  11, . 

The  first  step  in  applying  this  techniciue  is  to  partition  the  iidct  into  channels  atid  s<*etions 
of  channels  to  obtain  reasonably  uniform  values  of  r-  in  each  segment  and  uniform  cross 
sections.  The  procedure  is  illustrated  in  Section  VI  by  application  to  Masonboro  Inlet.  The 
next  step  is  to  select  a tentative  value  for  Manning’s  n.  The  values  of  /3,  F",  K,  and  r^  are 
then  evaluated.  At  time  equals  zero,  V and  II  may  both  be  taken  as  zero.  By  repeated  use 
of  equation  (23b)  to  define  equation  (41))  to  define  A and  integrating  equations  (25b) 
and  (26b)  for  one  time  step  values  of  V(t)  and  H,(t)  are  generated  from  prcscrilx'd  values 
of  tlie  open  side  tide  Il2(t). 

Calculations  for  the  tide  in  the  bay  are  rarely  satisfactory  for  the  first  few  tidal  cycles 
Ix'cause  the  proper  initial  conditions  are  unknown  and  the  assumption  that  II,  and  V are 
initially  zero  is  not  very  accurate.  However,  friction  causes  the  influence  of  the  initial 
conditions  to  decay  with  time,  and  the  computed  bay  tide  becomes  periodic  after  a few 
cycles.  When  the  agreement  between  computed  and  observed  bay  tide  becomes  acceptable, 
the  model  is  said  to  be  calibrated.  It  is  rare  to  obtain  satisfactory  agreement  with  the  first 
value  of  Manning’s  n selected. 

If  the  first  calculations  are  not  satisfactory  , a new  value  of  n is  selected  and  the 
procedure  is  repeated.  Several  values  of  n may  have  to  be  tried  before  satisfactory 
agreement  is  obtained. 
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APPENDIX  C 


DOCUMENTATION  OF  SECTION  PLOTTING  PROGRAM  SECPLT 

Program  SECPLT  is  used  to  plot  and  compute  areas  for  inlet  cross  st;ctions  from  digitized 
liydrographic  data.  Section  data  are  read  from  prepared  data  cards  and  the  program 
produces  a listing  of  input  data,  a listing  of  accumulated  cross-sectional  areas  at  specified 
width  increments,  and  a tape  of  plotter  data  for  CALCOMP  drum  plotting  of  the  section(s). 

Tlie  program  allows  variation  of  the  water  level  datum  for  area  calculations  by  allowing 
the  user  to  specify  the  number  of  the  data  point  which  should  be  used  as  the  vertical  datum. 
\\idths  tabulated  with  corresponding  areas  are  then  referenced  to  the  specified  point.  The 
final  tabulation  table  is  useful  in  determining  input  data  for  tlie  INLET  program. 

Description  of  input  and  output  follows;  the  program  is  coded  to  operate  on  the  WES 
Honeywell  437. 

1.  Input  Data. 

Card  Type  A. 

A 1— This  card  contains  the  number  of  sections  to  be  plotted,  the  number  of  parallel 
inlet  channels,  and  an  overall  scale  factor.  The  scale  factor  allows  the  plot  size  to  be  reduced 
by  the  factor  indicated  and  should  be  a number  not  exceeding  1.0.  FORMAT  (215,  F5.2). 

A2— This  record  consists  of  up  to  32  contiguous  alphameric  characters  which  will  be 
printed  out  to  identify  all  output.  FORMAT  (4A8). 

2.  Section  Data. 

A blank  card  is  entered  to  indicate  the  end  of  a data  set  for  a section;  an  additional  blank 
card  is  added  at  the  end  of  all  data  sets. 

Card  Type  B. 

Bl— The  first  data  card  for  each  section  contains  the  number  of  the  data  point  to  be 
used  as  the  datum  for  the  seetion,  the  horizontal  increments  at  which  areas  are  to  be 
eomputed  and  tabulated  (maximum  of  200  increments  per  section),  the  horizontal  scale  to 
be  ustid  in  plotting  (feet  per  inch),  and  the  vertieal  scale  to  be  used  in  plotting.  The  length  of 
the  y-axis  below  the  datum  is  fixed  at  6 inches  and  the  seale  should  be  chosen  accordingly. 
This  length  corresponds  to  a scale  factor  of  1.0  (card  Al)  and  will  be  reduced  if  a reduction 
factor  is  specified.  FORMAT  (15,  F6.0,  2F6.1). 

B2— This  record  is  a 24-character  label  and  is  annotated  on  the  associated  section  plot. 
FORMAT  (3A8). 

B3— Data  records  are  entered  for  each  data  point  in  this  section,  and  contains:  (a)  A 
point  number  which  provides  a notation  for  convenience  in  identifying  tlie  point.  A number 
must  be  entered  but  its  value  is  ignored  to  allow  for  introduction  of  additional  intermediate 
data  points  if  required,  (b)  horizontal  coordinate  of  the  point  in  feet,  and  (c)  elevation  of 
the  point  in  feet;  positive  values  are  below  the  zero  elevation  and  negative  values  are  above. 
FORMAT  (15,  F6.0,  F6.1). 
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3.  Output. 

As  each  x-section  set  is  processed,  several  diagnostic  notes  are  printed  to  help  in  locating 
data  errors  should  a run  prove  unsuccessful.  After  a group  of  up  to  six  sections  has  been 
processed,  the  output  routine  prints  out  the  output  tables  for  the  group.  These  tables 
include  a listing  of  tlie  input  data,  a listing  of  the  x-section  areas  as  functions  of  the  distance 
along  the  section  length,  and  a table  of  flow  grid  data.  The  flow  grid  data  may  be  useful  in 
establishing  the  INLET  program  data  for  determining  the  equivalent  prismatic  properties  of 
an  irregular  inlet. 

To  produce  the  grid  data,  tlie  program  divides  tlie  total  section  area  into  as  many  equal 
parts  as  the  number  of  parallel  channel  segments  specified.  Then,  tlie  cumulative  and 
incremental  widths  at  which  those  area  increments  occur  are  determined  and  printed. 

This  program  can  be  useful  but  must  be  used  carefully.  The  user  should  review  the  plots 
to  verify  that  they  are  reasonable  interpretations  of  the  data  points.  In  plotting  the  cross 
section,  the  program  uses  a cubic  spline  curve-fitting  technique  which  may  introduce 
unnatural  humps  in  the  section.  If  unusual  perturbations  occur,  they  can  usually  be 
eliminated  by  introducing  additional  data  points  in  the  given  region.  This  problem  is  mostly 
encountered  in  regions  where  the  slope  changes  rapidly.  When  a reasonable  plot  is  obtained, 
the  user  should  plot  the  flow  grid  over  tlie  inlet  plan  view  to  assure  that  the  grid  as 
determined  by  the  program  allows  for  a reasonable  flow  pattern.  If  the  pattern  is 
unacceptable,  it  will  be  necessary  to  adjust  the  grid  and  associated  data  before  using  the 
INLET  program.  It  may  be  found  that  the  x-sections  selected  cannot  produce  a reasonable 
flow  pattern  and  that  they  must  be  revised. 

The  manual  plotting,  planimetering,  and  data  calculations  for  the  x-sections  are  very 
tedious;  however,  initial  use  of  manual  techniques  may  be  helpful  in  developing  an 
appreciation  for  the  limitations  of  tlie  program. 


73 


SAMPLE  lOTUT  FOR  SECTION  PLOTTING  PROGRAM 


.JOB,WINTERGE  , RO  HM  803-G9R0-181  10 

.OPTIONS,,,,,? 

.SODA,V300,PRNT 

.PERIPH,T0404 

.LOAD 

(INSERT  OBJECT  PROGRAM  HERE) 

.DATA 

1 4 050 

MASONBORO  INLET  JULY-AUG  1966 


2 

20 

2000 

SECTION 

4— BAY 

END 

1 

-50 

-38 

2 

00 

-19 

3 

50 

00 

31 

200 

15 

32 

400 

18 

36 

800 

20 

38 

1000 

35 

4 

1300 

70 

5 

1360 

100 

55 

1600 

175 

6 

1800 

200 

65 

2000 

170 

7 

2250 

100 

8 

2350 

50 

9 

2410 

00 

10 

2470 

-38 

H52961HV11 
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SAMPLE  OUTPUT  FROM  SECTION  PLOTTING  PROGRAM 


MASONBORO  INLET  JULY-AUG  1966 

READ  SECTION  DATA 

X(l)  = -50.00  Y(l)  = -3.80 

X(N)  = 2470.00  Y(N)  = -3.80 

SECTION  DATA  READ  16  POINTS,  ORIGIN  = POINT  2 
WIDTH  INCREMENT  =20.0 
SCALE:  200.0  HORIZ 

5.0  VERT 

CALL  SPLINE 
SPLINE  COMPLETED 

INITIALIZATION  AND  AXES  COMPLETED 
K = 122,  KMAX  = 123 

MASONBORO  INLET  JULY-AUG  1966  L = 1 NUM  = 1 

KMAX  = 123  NUMWR  = 1 

OUTPUT  ROUTINE  REACHED 
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PAGE  1 MASONBORO  INLET  JULY-AUG  1966 
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(N 
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m 

tn 

NO 

00 

o 

(N 

m 

1-H 
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tH 

rH 

(N 

(N 

<N 

<N 
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SAIIPLE  OUTPUT 
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SAMPLE  OUTPUT 


INLET  FLOW 

GRID 

MASONBORO 

INLET  JULY-AUG  1966 

CUMULATIVE 

AREAS 

CHANNEL 

NUMBER 

SECTION 

1 

2 

3 

4 

1 

5869. 

11739. 

17608. 

23478 

CUMULATIVE 

WIDTHS 

CHANNEL 

NUMBER 

SECTION 

1 

2 

3 

4 

1 

1344. 

1689. 

1964. 

2450 

WIDTHS 

CHANNEL 

NUMBER 

SECTION 

1 

2 

3 

4 

1 

1344. 

345. 

276. 

486 
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NG  PRO 


MRSG.^'BO 


6BrjV4/MB[f  COPY 

OSTtS,  Sl‘<s  son  CABOS  V#V/I  I 

T*5K(Tnor37?uju,pwCEWCl  |TBT8)SI”S 

FTs. 

attac>‘('’LCT,pi.OTTT»  IO«OQT2u«n.^B«l  .cv«n 

LOAC(BLOT) 

LCO, 

BE«I'-Of  T4PE5) 

C0P»S6B( Tap? 3 •OUTPUT) 

pboc-Ba"  secplT( input, Output, Tapes* I^PuT,TAPE6«nUTPUT I T apes • punch* 
1TAPE3) 

C SECTION  plotting  PBOGRap 

C OaPS  VEPSION--OCT  71 

c 

c 

OI-FNs’PN  IbUFPdOOO) 

0I*<ENsI0N  7(100),  V(ion),  A(IOO),  TlTLl(<J)t 
2TtTL?T5)»  '-RTARr(2,b«2no)  • D A T A ( 2 , fc  , 1 00  ) 

''I“ENstON  04(10),  •■(fe*10)«  OAaOO) 

OATa  aXSLBL/ShOATU"/,  TP01,Iwr1/5,6/ 

OATa  vIBL/UhELEV/ 

OATA  jolTPE/  3/ 

PE.jNn  IPLTPE 

CALL  PLOTS  (leUAE,  100n,IPLTPE) 

CALL  PI nT(o,*0,,«3) 

C BEAn  C4Rr^  TYPE  A1  containing 

C ‘-PLTS  (IS)-  nUkheb  of  sections  TO  BE  PLOTTED 

c ‘cnnl  (ts)-  number  of  Parallel  inlet  Channels 

c SCALE  (Fs,2)-  overall  Scale  factor 

<*  BEAU  (TRDl,200  t)  nPlTS,  nChNL,  SCALE 

IKE  OF  ( IBOl  ) )<»R,R 

e IF  (nplts)  rr,  rr,  r 
R NC-l  * nCmnl  - 1 

c 

c N » Pt  pTTeo  section  index 
c 

Nip 

CALL  factor  (SCALE) 
vPO  I 10, 

C RFao  CaRo  type  a?  containing 

C TITLl  (UAR)»  32  CmaRACTFB  OuTPut  I OFNT I F IC AT  I ON 
read  (tROl,1001)  TITLl 
.RITE  (IrR1,1028)  TITLl 
«PP  I ?, 

c 

C NUN  * absolute  index 
c 

NL“  « 0 

»“A«  a 0, 

OU  OS  L « ' • nplts,  N 
•'L  2 T « t , R 
00  1 J a 1 , 1 00 
■'ATa(i  , I,  J)  « 0, 

1 f'AT4(?,l.Jl*0, 
no  2 J 1 1,  200 
.WT4R»( 1 , 1 , J)  ■ 0, 

2 -hTaRv(?,I,J)  « P, 

C 

c clear  array  after  each  sIX  SFcTIOnS 

c 

NM4X  I 0 


82 


8 0 

NU*^V.«  ■ 0 


besi:avaiiable_cop)( 


c 

C NU’^kO  ■ xHITE  SET  INOEx 

C 

s .BITE  fi.ptitoin 

C ME*n  c*Bn  type  Bi  containing 
C (15)-  nu“BEb  {jf  data  point  to  «e  used 

C AS  TH,r  oath.  AUP  T^E  SECHON 

C 0*  (^■^,0)-  mCjRIZOnTAU  InCBEmENTS  At  .HlfH 
C ABEAS  ABE  TO  BE  Cn'"BUTED  anO  TaBULATEO 

C SCALE*  (Pfe.n-  hobiZONTau  scale  to  re  USEO  In 
C Plotting  (AEET  PEB  Inc**) 

C 5CALV  (P6.11-  VERTICAL  SC*LE  Tp  BE  iJSEO  IN 
C Plotting  (PEET  per  iNCR) 

BEAO  (TB01.200R)  NZ£»«  PX»  SCALXi  SCALV 
C BEAO  CA«n  TYPE  82  containing 

C TITL2  (JAR)-  2u  C«aBaC1FB  LABEL  ANNOTATED  OH  TMf 

c associated  section  Plot. 

BEAO  (IBDl.lOon  TnL2 
NPTS  c I 

C READ  CABO  Typt  BJ  CONTAINING 

c j (IS)-  POINT  number 

C X(nptS)  (Ee.O)-  horizontal  C0-0Rd1NaTE  OE  POINT  IN  FEET 
C y(NPTS)  (Fb.D-  elevation  of  point  tn  feet 
C mEao  One  card  for  each  data  point 
C LAST  CaBo  CP  OaTa  sft  EOR  a SFCTIOn  is  blank 
c Of.F  AOntTIONAL  blank  CARo  at  ENp  OF  all  OATA  SETs. 

10  REAP  (TPOl,200«)  J.  X(NPTS).  V(NPTS) 

IF  (J)  IS.  IS.  12 
12  NPTS  > NPTS  ♦ I 
GO  TO  IP 

IS  NPTS  a nPTS  - 1 

-RITE  (I.«I*1020)  X(l),  V(l),  X(NPTS),  V(nPTS) 

.RITE  (Iwwi.lOOB)  NPTS.  nZeR.  OXt  SCALX.  SCALY 
kChk  a 0 

DO  25  J » 2t  nPTS 
IF  (XrJ)  - X(J  - J))  2u.  2<H  25 
2u  -ritp  (I-ri.1030)  J-li  x(J-l).  X(J) 
kCmk  a KCmk  ♦ 1 
2S  continue 
VOO  » n. 

xTBn  a X(NZE«)  - X(l1 
.RITE  (I-R1.1012) 

IE  (kcmkI  2T»  27.  28 
C 

C OETeR-TnE  spline  constants 

c 

2T  call  SPLINP  (X.  V,  nPTS,  A) 

.RITE  fl.Hl.IOlJI 

c 

c SPLINP  COMPLETED 

C 

2*  XT  ■ x(l) 

0X5  « P,5»0X 

area  a 0, 

SCAL2X  « SCALX  ♦ SCALX 
SCXINy  a 1. /SCALX 
• SCVlNy  a -1  ./SCALY 

XP  a x(I)»SCXINV 
VO  a y(l)*SCVINV 
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v*0  • Y(nZt«)*SCVINV 
vtpn  ■ r(NZER) 

VV1  ■ v»n  ■ 6,0 
XJ  ■ X(NZEB)*SCXINV 
XXt  ■ .2»IHX(0,S*9CXlwy*XTBK'  « 1,8) 

XSOTE  • SC*1.*»*X1 

XXI  ■ y|  * XXI 

XAO  > ?,  « SCXlNV»X(1) 

*L  « scxiNv»(x(NPTS)  • xtn)  ♦ «,o 

C*LL  PLOT(xPO,0,,.X) 

C41.L  PLOT  (X*0t  vPn,.J) 

c 

C PLOT  X AXIS 

C 

CALL  AXlS13(XXl«VXfliXXSLBL,Si9,aiiXLT>t •OiXMOTEtSCALaXi 

X2,0,«i ,0) 

C 

C PLOT  nATX  POINTS 

C 

00  SO  I « ti  nPTS 
OATaI 1 ,NUMW«T J 1 1 ) » x(l) 

OaTaIJ.NUMkHTI*!)  • V(T) 

IS'JB  ■ nPTS  • I ♦ 1 
XP  s SCXINV»X( ISUB) 

VP  « SCVINV*V{ ISUH) 

TP  (XP  . mo,)  31  i S?t  32 
SI  IP  (VP  « 30,)  33t  32*  32 
52  KNlTfc  (IVH1  * 1032)  ISUb,  X(ISUB)*  VdSUB) 

<CNx  z kCmk  ♦ I 

GO  TO  30 

33  CALL  SVHBOL  {XP*VP.0,12*S,0,*«I) 

30  CONUk'IiE 

C 

C BLOT  titles 

C 

CALL  SX^BOL  (0,*-T,.0,28.TITLlt0.iS2) 

CALL  SX-BOL  (0, .-7,5.0, 28, TITL2,0,,2«) 

SCXL2V  ■ SCXLV  ♦ SCXLT 
VnOTE  ■ •6,*SC*LX 
C 

C PLOT  y AXIS 

C 

CALL  AXISIS  (XI, VYt.VLBL.U. 0,21,8,01-1,1, VNOTE, 
1SCXl2v,2,0,-1,0) 

IP  (XCMK)  35,  35,  72 

C 

C “OVt  TO  PIPST  POINT 

C 

35  CXLL  PLOT  (XP,XP,3) 

NIJ-  * non  ♦ 1 
NU-kP  ■ NUMWP  V t 
X ■ 0 

"RITE  (I"R1,101«) 

C 

C AXES  COMPLETED 

C 

JPT  ■ 1 

50  XT  z XT  ♦ OX 

C 

C INTtPpOLXTE  POR  V(X) 

C 


BEST  AVAILABLE  COPY 

C*LL  SPLlNT(XTtYOiVl •X,Vi*,NPT8t JPT) 

C 

C compute  IP  V BELOW  *KIS 

C 

TP  (Vfl  • YTHN)  60f  bOt  55 
55  AREA  ■ AREA  * DX5«(Y0U  ♦ Yn  ■ YTRN) 

« ■ < ♦ I 

WRTARyT t iMIWWNiK)  ■ XT  ■ X(NZER) 
wrt*bv(2«mi^»R'B)  ■ area 
YOO  ■ VO  • YTRN 
60  XP  ■ SCXl^'YAXT 

YP  ■ SCYInvPYO 

C 

C PLOT  SECTION 

C 

CALL  PLOT  <XP.YP,2J 
IE  (XT  • X(NPT8)  ♦ OX)  5O1  50,  70 
70  wRTxRyf 2#M)«wRiK*l ) ■ AREA  ♦ 0X5*700 

wMTaRv  ( 1 I NLiMxRtXA  I ) ■ wrT  APV  ( I iRURwRi  X)  A OX 
nATMjMhW)  M »RTARY(2*NI)MrtH,KA)  )/NCHNL 
W(NU*'i.BtNCHNL)  ■ wRT  ARY(  1 tNUN«»#KA  I ) 

C 

C CLOSE  ON  last  data  point 

C 

XP  ■ 8CXINV»X(NPTST 
YP  ■ SCyInv»Y(NPTS) 

CALL  PLOT  fXP,YP,2) 

72  IE  (N-AX.LT.npTS)  NmAX  • NPTS 
lA  (KPAX  • K • 1)  71*  7U*  7« 

Yl  WMAX  ■ K * 1 

-MITE  (I-R1*I0lb)  K,  KmaX 
7a  N ■ N A t 

lA  (NV)«  . NPLTS)  75.  'AO.  50 
T5  if  (N  . 2)  eo.  85.  65 
60  XPO  S -XAO 

XMAX  ■ XP 
YPO  ■ 2U, 

GO  TO  5 

65  IE  (Xmax  • XP)  86.  eS.  66 
86  X-AX  s XP 

86  XPO  • XMAX  A 5.0 

VPO  ■ 10. 

N » 0 

IE  (N(imhR  . 6)  5.  60*  60 

60  WHITE  (IWRI.IOIO)  TlTLl .L»NUM,KMAX,NUmwB 

WRITE  (IwRl.1015) 

IPG  M NMAX/50  A 1 
00  6?  .ti  ■ 1.  IPG 
C 

C output  routine 

C 

.mite  (IwMt.1017)  J3.TITL1 . ( J. J«L.NU“) 
white  (IwHl.lOlH) 

T-Ax  ■ 50*Ji 
T“IN  « I“*X  - «9 
IF  (Js.EO.IPG)  IMAX  B NhAx 
00  62  M I-IN.  IMAX 

«2  white  ( I-Hl  . 1016) (0ATA( 1 . J1 . J2) .DATA(2. J1 . J2) . JI»I.NUMWH) 

IPG  « Kmax/50  a 1 

00  I65  JO  • 1.  IPG 

.RITE  (I-H1,1005)  JU.T TTLl . ( J. JbL.NUH) 
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-HITE  (I-Hl*t004) 

IHAX  s bO*JU 

IMJn  g JMAX  • U4 

IE  (Jo. to, IPG)  IMAX  a XMAX 
DO  1<>S  J?  a I-IN.  IMAX 
IF  - n 19!>i1P5,91 

9l  no  9«  J3  a It  NCHl 

DO  94  Jt  * I*  ^UM-P 

OOA  a FLOAT(J3)anAf Jt ) 

IF  («pTAHYf2iJUJ2)  • no*)  93.  93.  9« 

93  IF  (-RTAHY f 2. J1 . J2* 1 ) « OOA)  94.  94.  193 
193  OX  B wHTAHy(l.Jt,J2An  . -RT ARY ( 1 . J I . J? ) 

DEL  a nx»(OOA  . ..HTARY(2.  Jl  . J2)  )/(-RTARYf2.Jl  . J2*n  • 

1-RTAHy(2. J1 . J2)  ) 

-(Ji.«J3)  a "HTAWY(  I ,J1  ,J2)  ♦ del 
9«  CONTI. lYlf 

mRITE  (l-Rt.lOOb)  (KRTARY(l.Ji.J2).-RTARV(2.Jl.J2)tJlBl.NUMHR) 
195  continue 

c 

IF  (NCHNL  - 1)  95.  95.  105 

105  WRITE  fI-R1.102n  TITLl.IJli  Jial.  NCHNL) 

DO  lip  Jl  a 1.  NCHNU 
UO  DAA(Jl)  a 0. 

DO  l3n  Jl  a 1.  nUmwH 
DO  115  J2  a 1.  NCHNL 
115  OAACJ?)  a J2*UA(Jl) 

ISO  white  (I-H1.1025)  Jl*(nAA(j2).  J2  a I.  NCHNL) 

WHITE  fI»Hl.1022)  ( J 1 • J 1 ■ 1 . NCHNL ) 

DO  135  Jl  a 1.  Nijw-H 

>35  WHITE  fI-Rl.l025)  Jl.  fwCJi.Ja),  J2  • 1.  NCHNL) 

WRITE  (IWRI.1024)  (Jl.Jl  a l.  NCHNL) 
no  145  Jl  a 1.  NUMWR 
OAA(l)  a -(Jt.l) 

00  140  J2  a 2.  NCHNU 
lUO  DAA(J;)  a wCJl.J?)  . w(Jl.J2«i) 

>45  white  fIwRl.1025)  Jl.  tO**(J2).  J2  a i,  NCHNL) 

95  continue 

00  TO  4 

c 

99  XHAx  a XMAX  ♦ 5, 

CALL  PLOT(XmAX.0,|999) 

RE-INo  IPLTPE 
STOP 
C 

lOOl  F0RHAT(3Alft.A2) 

>004  FORMAT  (6(7X.  iHW.  9X.  IHa.  2X)/) 

>005  format  t5MlPAOE.I3.40X.4Ae//27H  SECTION  NUHBEHSi  AREAS-KSF/ 
t6(12X.  13.  5X)) 

1006  format  (bCpPFlO.O.  mIPFIO,]!)) 

1006  format  (/IflM  SECTION  DATA  READ  14.  9”  POINTS.  . 

HRHORlGIN  a POINT.  I4/16H  WIDTH  INCREMENT  ■.  F?,!/ 

»7H  SCaLEI.  F»,1.  6M  HORIZ/  7X.  F6,l.  SH  VE«T) 
loto  format  (/IX.  4Ae.  5X.  3HL  a.  13.  SH  HUM  a,  !«/ 

17H  KM*x  a.  14.  lOH  NUMWR  a.  U/) 

1011  format  (18H  READ  SFCTION  OaTA  ) 

1012  format  (12H1CALL  SPLINE) 

1013  format  (ITH  SPLINE  CQhplETED) 

1014  format  (34h  Initialization  and  axes  completed; 

1015  format  (21H  OUTPUT  ROUTINE  REACHED/) 

1016  format  (4h  K a.  I4.  8H.  Kmax  a.  I4) 

lOlT  format  (5h1PACE.I3i40X,4A6//15H  INPUT  dATA—FT, 


f 
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BEST  AVAIWBl!  COPY 


lote 

1019 

1020 

1021 


1022 

102<t 

1025 

1020 

1030 


1032 


C 

2001 

200tt 


C 

c 

c 

c 

c 


c 

c 

c 

c 

c 


30 


C 

c 

c 


60 


best  AVAIUBIE  COPY 

FOMMAT  C0(TXf  IHXt  IHy,  2X)/) 

format  (fcCFlO.Oi  Ft0,2)) 

F0»mat(7m  x(1)  ■«  F6,2,  9H  Yd)  at  FS,2/ 

ITM  X(M  »*  F0,2.  9H  V(N)  a,  Fe,2) 
format  dbMllNLiT  FLOm  CRIOi  lOX*  aA0// 

M7H  Cumulative  areas/ 

Aaoxi  iSM  CHAMNtL  mumbeb/  an  section*  10*  10112) 
format  t/lRM  cumulative  mioTMS/  MOX* 

HaMCMANNFL  NUMBER/  6H  SECTION,  16,  101)2) 

format  (/7h  »IDTmS/uO*,  IrmCHannEl  NUMBER/ 
section,  Ib*  10112) 
format  (15.  11E12,0) 
format  (1M).  «A8/) 

format  {/uom  •*«*«a*«**a**a**««*«// 

2 t6M  ERROR*  ASSISSCA*  (U,  (H,,  Fb.O* 

3 28H  EXCEEDS  SUCCEEDING  VALUE  OE*  Fb.O// 

tt  uoH  aa««**«aaa«a«*»*«a*«/) 

format  (/ROM  »****a»**a***»**»»**// 

2bM  POjNT*  I«*  2FB,?*  2ih  OUT  OF  PLOTTER  RANGE// 

3aOH  aaaa««aa«*Aa*aaa««a*/) 


format  (2I5.F5.2) 
format  (15*  Fb.O*  2Fb.i) 

END 

SUBROUTINE  SPLINE(X*V*N,Y1) 


REAL  L 

N a Number  of  points*  eouations 

VJ  a calculated  first  OE«IVATIVtS  OF  V 

DIMENSION  X(N),  V(N)*  V1(N)*  A(IOO)*  BdOO) 

N1  a N • 1 

L ■ 1,/AHS(X(2)  - X())) 

01  a 3,A(Y(2)  • V(n)*L*l. 
n(l)  a L 

A(l)  a B(l)  * B(l) 

Vl(l)  a Ol 

A a vector  of  principal  Diagonal  elements 
R a vector  of  1ST  OFF  DIAGONALS 
STORE  CONSTANTS  In  VI 

00  50  I ■ 2*  Nl 

L ■ 1 ./ABS(X(  !♦! ) - X(  I) ) 

02  « J,»(V(IA1)  • V(I))*L»L 
n(I)  a L 

A(I)  a 2,*(B(I)  * R(I«1)) 
ri(I)  a 01  ♦ 02 

01  a 02 

A(N)  a 8(N1)  ♦ B(NI) 

V1(N)  a 01 

equations  set*  symmetric  and  TRIDIACOnaL*  nor  triangulate 

00  bO  I a 2,  N 
F a B(T-1)/A(I-1) 

A(I)  a A(I)  <*  B(Im1)*F 
Yl(I)  a Yl(I)  • Y1(Im1)*F 
Y1(N)  a Y1(N)/A(N) 
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BES1'AVAIIA31E  COPY 

c 

C SOLVE  BY  HACK  substitution 

c 

do  70  I ■ If  Nl 
J ■ N . I 

To  Y1(J)  ■ (Yt(J)  • B(J)*Y1(Jy1))/*(JI 

return 

END 

SUHHOuTiNt  SPLINT  f XOivOfVlOf XfYf Y1 |N, J) 

C 

dimension  V{n)i  Y(N)f  Y1(N5 
C 

C J indicates  span  of  last  call  to  routine 

C if  X ,LT,  X(J)i  start  at  beginning 

c 

IF  (Xo  • X(J))  lOi  aoi  20 
10  J • I 

20  00  2S  I ■ Ji  N 

IF  (Xfl  . X(I))  <tO,  60|  25 
2b  CONTINUE 

c 

c xo  ,GT.  X(N)t  extrapolate  FROM  X(N) 

C 

J ■ N 

YO  ■ y(M  ♦ IXO  • X(N)T»YlfN) 

YIO  ■ Yl (N) 

return 

«0  IF  (I  - 1)  <45i  Rbf  50 
C 

c XO  .LT.  x(i)f  extrapolate  from  x(n 

c 

tt5  YIO  ■ Yl  (n 

YO  ■ Y(l)  ♦ (XO  « X(1))*Y10 

return 

c 

C interpolate  on  span  j 

c 

50  J ■ I - 1 

TL  • 1 ./(X( J*1 ) • X( J) ) 

xoo  ■ xo  • x(j] 

C • Y 1 { J ) 

0  ■ Y(J) 

8 ■ (YfJYll  • T(J))*TL*TL 
A « Tl*(TL»(Y1  (J)  4 Ylfjtl))  . B . B) 

8 ■ J,»B  - TL*(Y1(J)  4 Yl(J)  4 Yl(j4n) 
YO  ■ A*x00»*3  4 R*X00*x00  4 C»XOO  4 0 
YlO  ■ 3,*A*X0C*X00  4 2.«B*X00  4 C 

return 

c 

C xo  IS  AT  A DATA  POINT 

C 

60  YO  ■ Yd) 

YlO  ■ Yl(I) 

J«I 

return 

END 


1 tt  050 

masonboro  Inlet  jult^aug  irbs 

2 20  2000  50 

Section  a>«nAY  end 
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1 

•5fl 

-38 

2 

00 

•19 

i 

50 

00 

J1 

200 

15 

32 

uoo 

18 

38 

800 

20 

38 

1000 

35 

U 

1300 

70 

5 

1360 

100 

5S 

1600 

175 

6 

1800 

200 

65 

2000 

170 

7 

2250 

100 

8 

2350 

50 

9 

2010 

10 

*4 


BEST  AVAILABIE  COPY 

10  2«T0  -58 
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APPENDIX  D 


DOCUMENTATION  OF  PROGRAM  INLET 

INLET  is  a digital  computer  program  whicli  performs  various  numerical  computations 
and  equation  solutions  required  in  the  lumped  parameter  analysis  of  an  inlet  bay  system. 
The  program  operation  can  be  divided  into  three  principle  functions:  (a)  Analysis  of  inlet 
grid  data  for  nonprismatic  inlets  to  determine  equivalent  prismatic  inlet  properties  according 
to  the  method  described  in  section  III.  This  function  may  be  deleted  if  not  applicable;  (b) 
numerieal  solution  of  tlie  simultaneous  nonlinear  differential  equation  (25)  based  on  input 
data  and  results  of  function  (a);  and  (c)  generation  of  a solution  plot  for  function  (b).  This 
feature  may  also  be  deleted  if  only  a tabular  solution  is  desired. 

1.  Input  Data. 

Input  for  INLET  should  be  on  punchcards  in  the  formats  shown  for  A and  B. 
Descriptions  of  individual  card  requirements  follow. 

Curd  Type  A. 

Al— This  card  indicates  whetlier  solution  plots  are  desired  at  any  (joint  in  the  program 
execution.  Al  in  cc2  of  this  card  indicates  that  plots  arc  desired;  if  no  ()lots  are  desired,  the 
card  should  be  left  blank. 

A2— This  card  defines  grid  and  system  data  for  the  chosen  flow  grid,  inlet,  bay,  and 
tidal  conditions  for  the  area.  The  first  value  indicates  the  number  of  flow  cliannels  in  the 
assumed  flow  grid;  the  second  value  similarly  indicates  the  number  of  x-sections  whicli  have 
been  selected  across  tlie  inlet.  If  the  number  of  x-sections  is  entered  as  zero,  the  inlet  is 
assumed  as  (irismatic  and  all  parameters  are  entered  on  A4.  The  next  two  values  arc  the  tidal 
[leriod  in  hours  and  the  tidal  semirangc  in  feet.  Manning’s  n for  the  inlet  is  specified  next: 
however,  if  it  is  desired  to  s()ecify  distinct  values  of  n for  each  channel  and  x-section,  zero 
should  be  entered  here;  the  distinct  values  will  tlien  be  read  later.  The  MTL  surface  area  of 
the  bay  in  square  feet  is  input  in  scientific  notation  with  the  mantissa  entered  to  the  left  of 
the  “E”  and  the  power  of  10  to  tlie  right.  The  dimensionless  surface  area  variation 
parameter  ^ follows  and  is  determined  as  previously  described.  The  final  value  for  the  card 
is  the  lieacli  slope  f defined  as  horizontal  to  vertical  differences. 

A3— These  two  cards  are  title  cards  and  each  may  contain  up  to  32  characters  of 
numerical  or  alphabetical  information  to  identify  the  resulting  output.  These  titles  will  also 
appear  on  the  solution  plots  if  any  are  generated. 

A4— This  card  is  omitted  if  equivalent  prismatic  properties  are  to  be  eomputed  from 
data  on  card  type  B.  This  card  should  be  used  only  if  the  number  of  sections  on  card  A2  is 
specified  as  zero,  otherwise  the  B cards  should  be  used  in  its  place.  For  an  inlet  of  known 
prismatic  properties,  the  data  values  on  this  card  include  the  MTL  cross-sectional  inlet  area 
(s<(uare  feet),  inlet  width  (feet),  inlet  depth  or  hydraulic  radius  (feet),  and  the  inlet  length 
(feet). 
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A5— This  card  sets  up  the  equation  solution  parameters.  The  first  value  is  the  time  at 
whieh  the  calculation  is  initialized  and  the  second  value  is  the  time  at  which  output  is  to 
begin.  If  it  is  desired  to  start  output  at  time  zero,  the  beginning  time  should  usually  be  some 
negative  multiple  of  the  tidal  cycle  in  hours,  thus  allowing  a transient  solution  stage  during 
which  the  effects  of  inexact  initial  conditions  may  be  dissipated.  The  third  data  value  is  the 
time  in  hours  at  which  the  solution  is  to  be  terminated. 

The  next  two  values  are  the  initial  conditions  for  bay  surface  elevation  and  inlet  velocity. 
For  sinusoidal  ocean  tides,  tbe  initial  bay  level  (negative)  may  generally  be  taken  as 
one-fourth  of  the  tidal  semirange  and  the  initial  velocity  will  likely  approximate  3 feet  per 
second.  White  the  initial  values  chosen  should  not  usually  be  critical,  tlie  closer  they 
approximate  the  steady-state  conditions,  tlie  shorter  may  be  the  transient  stage  of  the 
solution.  The  next  value  specifies  the  time  step  in  minutes  to  be  used  in  the  numerical 
solution  method.  A 5-minute  time  step  has  been  used  successfully  for  sinusoidal  tides  but 
may  be  varied  if  desired.  All  time  steps  need  not  be  output  as  tlie  seventh  data  value 
specifies  which  values  are  to  be  printed.  For  example,  if  a 5-minute  time  step  is  entered  and 
a 12  is  entered  for  the  number  of  steps  between  outjiut,  output  will  be  printed  every  60 
minutes.  Critical  values  of  all  functions  will  always  be  printed. 

The  final  value  is  a plot  indieator  and  is  set  eijual  to  one  if  a solution  plot  is  desired.  If 
left  blank,  no  plot  is  generated  and  tins  card  terminates  the  data  set. 

A6— This  card  describes  the  data  which  establish  the  limits  and  scales  for  plotting  the 
com[)uted  solution  if  such  a plot  is  desired.  The  first  two  values  define  the  starting  and 
ending  times  for  the  plot  and  should  not  exceed  the  range  of  the  tabulated  output  specified 
on  eard  A5. 

The  time-axis  scale  (hours  per  inch)  is  specified  next  and  should  be  chosen  to  result  in  a 
plot  length  of  at  least  10  inches  to  complement  the  legend  and  title  format.  For  the  vertical 
axis,  the  minimum  value  of  tidal  heights  (feet)  or  velocities  (foot  per  second),  overall  height 
of  plot  (inches),  scale  of  tidal  heights  (foot  per  ineh),  and  velocities  (foot  per  second), 
minimum  value  of  flows  (thousand  cubic  feet  per  second),  and  scale  of  flow  (thousand  cubic 
feet  |)er  second  per  inch)  are  specified. 

The  final  data  value  on  the  card  allows  adjustment  of  the  total  plot  size  by  tbe  value 
selected.  For  instance,  if  a single  sheet  plot  size  is  desired,  the  above  data  should  be  selected 
to  produee  a plot  of  approximately  5 by  10  inches.  Specifying  a size  factor  of  0.75  will  then 
produce  a plot  3.75  by  7.5  inches  which  will  fit  on  a standard  sheet  size  with  axes  and 
annotations.  The  size  factor  could  also  be  used  to  increase  the  plot  size  if  desired. 

2.  Irregular  Inlet  Grid  Data. 

Curd  Type  B. 

Bl  — Label  card  to  indicate  that  data  following  are  cross-sectional  areas. 

B2— Area  eards;  each  card  contains  tlie  cross-sectional  area  of  each  assumed  flow 
channel.  One  record  should  be  used  for  each  cross  section. 
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B3— Label  card  for  widths. 

B4— Width  cards;  each  card  contains  the  widtli  of  eacli  assumed  flow  channel  for  the 
corresponding  section. 

B5— Label  card  for  lengths. 

B6— Lengtii  cards;  each  card  contains  tlie  lengtli  between  cross  sections  along  the 
boundaries  of  the  assumed  flow  channels.  The  number  of  lengths  per  card  is  one  greater 
than  the  number  of  cliannels  and  tlie  number  of  cards  is  one  less  tlian  tlie  number  of 
sections. 

B7— Label  card  for  Manning’s  n. 

B8— Values  of  Manning’s  n for  each  channel  of  each  cross  section.  A value  sliould 
correspond  to  each  of  the  B2  and  B4  data  values. 

Cards  B7  and  B8  must  be  omitted  if  a value  of  Manning’s  n otlier  tlian  zero  is  specified 
on  card  A2. 

Cards  A2  through  B8  may  be  repeated  if  several  problems  are  to  be  solved  in  a single  run. 
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(Utmirtfr«4uftiomofthUformUmithoHie4c/.frdfirf4informa1iomttiiafdi>t.>  tO  COLUMN  KFY  PUNCH  TRANSCRIPT  LAYOUT  SHIET 


r AD-A0b2  600 
UNCLASSIFIED 


COASTAL  CN6INEERIN6  RfSCA'^CH  CENTER  FORT  BELVOIR  VA  F/6  8/8 

comparison  of  numerical  a 40  PHYSICAL  HYDRAULIC  MODELS*  MASONBOR— ETC (U) 
JUN  77  C J HUVAL*  « L HINTER6ERST 

CERC-8ITI-8-APF-8  NL 


END 


DATE 

FILMED 
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SAMPLE  INPUT  FOR  INLET  PROGPJUI 


.JOB.WINTERGE  RO  HM  803-G9R0-182  10  H52962HV1 1 

. OPTIONS,,,,, F 

.SODA, 0300, COMP 

,PERIPH,T0404 

. DFAREA,00 ,AM,ANY , 2000 , 48 

.LOAD 

(PROGRAM  OBJECT  DECK  INSERTED  HERE) 

.DATA 

1 

4 5 124166667  1900  00270  1866&8  03500  7500 

MASONBORO  INLET  1969 
MEAN  SEA  LEVEL  CONDITION 
AREAS 


8350 

8290 

8460 

8280 

4990 

4240 

6710 

3950 

3090 

6600 

4770 

4820 

4010 

4400 

3370 

2320 

4550 

3220 

4510 

3280 

WIDTHS 

1080 

870 

1050 

725 

905 

995 

475 

330 

550 

190 

150 

360 

500 

130 

180 

345 

LENGTHS 

960 

1100 

1200 

1440 

1500 

650 

1070 

1290 

1510 

1640 

190 

540 

630 

700 

870 

350 

525 

540 

550 

500 

-1241666 

000000 

1250000 

-0500  3000 

5 6 

1 

000000 

1242050 

125000 

-4000  4800 

1667-60000 

25000 

15 
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SAMPLE  OUTPUT  FOR  INLET  PROGRAM 


EQUIVALENT  COEFFICIENT  OF  REPLETION 
FOR  HON-PRISMATIC  INLET 

MASONBORO  INLET  1969 
MEAN  SEA  LEVEL  CONDITION 


AREAS 


8350. 

8290. 

8460. 

8280 

4990. 

4240. 

6710. 

3950 

3090. 

6600. 

4770. 

4820 

4010. 

4400. 

3370. 

2320 

4550. 

3220. 

4510. 

3280 

WIDTHS 

1080. 

870. 

1050. 

725 

905. 

995. 

475. 

330 

750. 

255. 

250. 

625 

550. 

190. 

150. 

360 

500. 

130. 

180. 

345 

LENGTHS 

0. 

0. 

0. 

0 

1030. 

1150. 

1320. 

720 

860. 

1180. 

1400. 

755 

365. 

585. 

665. 

350 

437. 

532. 

545. 

275 

0. 

0. 

0. 

0 

TIDAL  PERIOD  = 12.417  HR 
TIDAL  SEMIRANGE  = 1.900  FT 

BAY  AREA  (MSL)  = 1.86600E+08  SQ  FT 
MANNING’ 'S  N = 0.0270 

COEF  OF  REPL  = 1.681 

INLET  AREA  = 1.41000E+04  SQ  FT 
INLET  LENGTH  = 3.04250E+03  FT 
INLET  DEPTH  = 11.591  FT 
INLET  WIDTH  = 1250.  FT 
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SAMPLE  OUTPUT  FOR  INLET  PROGRAM 


INLtT  ANALYSIS 


MASON80RO  INLET  1969 
MEAN  SEA-  LEVEL  CONDITION 


TIDAL  SEMIRANGE  = 1.900  FT 

MFAN  HAY  SURFATP  ARPA  r 9 n FT 

BAY  SIDE  SLOPE  PARAMETER  BETA  s 0,350 


inlet  properties: 

Y-SPrTlON  ARFA  RFI  nu  MTi  s l Amnp^QA  <;q  rT 

INLET  WIDTH  AT  MTL  = 1250,  FT 

INLET  hEACH  SL-AP-F  = iSU 

MTL  DEPTH  = 11,591  FT 

LNLEJ-LENGTH  = 3042.  FI 

MANNING'*S  N = 0.0270 


T IMF 

I NFl.  QU 

HI 

VF| 

n 

HRS 

FT 

KCFS 

F T 

FPS 

KCFS 

u.oo 

2,82 

39,04 

n . “jn 

HKWnVBH 

44 . 22 

1.00 

0.92 

0.41 

3.20 

47.83 

1.17 

1 . i;6 

n.55 

48.54 

1.50 

1.31 

0.00 

0.83 

3.16 

49. 05» 

p . iin 

1 . fci 

11 . n II 

1 . ?? 

2.0.3 

46,07 

2.50 

i.ei 

0.00 

1.56 

2.47 

40.50 

3 . nn 

1 . vn 

1 .82 

1.69 

28.1  a 

3.08 

HBKffl 

0.00 

1.85 

1.52 

25,38 

.3  . >3(1 

1 .95 

n..3a 

6 , 36 

3,58 

IBQfl 

0.00 

1.95» 

0.06 

0.96 

4 . on 

1 .87 

-1.58 

• 9 /i  . i A 

4,50 

ID^B 

0.00 

1.65 

-2.03 

-32,87 

1.39 

-2.41 

-37.93 

5 ,50 

0.67 

0 .00 

1.07 

-2.80 

-42,  75 

5.,.0  0_ 

Q-.-20_ 

0.00 

0 .71 

-3.11 

-45.70 

6.42 

-.20 

0.00 

0.37 

-3.27 

-46.50* 

-.28 

II . 3 r 

-3.29 

-46 . 48 

6.83 

-.59 

0.00 

0.01 

-3.33# 

-45,80 

-.7  4 

HiVilliTSH 

-.14 

- .1  . .1 .3 

-45 . 1 1 

7.50 

-1 .16 

0.00 

-.58 

-3.19 

-41,71 

a.un 

-1  .50 

-1.03 

-2.87 

-36 . 33 

8.50 

-1.74 

0 ,00 

-1.44 

-2,35 

-28,91 

HHHKJviTil 

-1  .88 

HHiWiTiH 

-1.77 

_ 1 _ S 

• 1 Q . M .T 

9.33 

-1.90* 

0,00 

-1.92 

-.87 

-10.48 

9.5n 

-1  .89 

-1.96 

-.43 

-5.20 

9.67 

-1.87 

0.00 

-1.97* 

0.11 

1.31 

10.00 

-1.79 

0.00 

-1.89 

1 . OH 

12^,96 

10.50 

-1.57 

0,00 

-1,68 

1.30 

15.92 

11.0  Q_ 

-1.25 

0.00 

-1.44 

1.67 

20.96 

11.50 

-.85 

0.00 

-1.13 

2.13 

27,55 

-.40 

-.77 

2.54 

34.01 

12.50 

0.08 

0,00 

-.37 

2.87 

39,96 

AVG  BAY 

LEVEL  « 

0.02  FT 

• CRITICAL  POINT 

VALUE 
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SAMPLE  OUTPUT  FOR  INLET  PROGRAM 


***  TABULATION  COMPLETED  *** 


ST  20 
ST  30 

PLOT  X AXIS 
PLOT  RIGHT  AXIS 
PLOT  LEFT  AXIS 
BORDER  COMPLETE 
50  I = 1 
80  I = 1 
80  I = 2 
80  I = 3 
ST  80 
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INLET  l>:U)Gll\:i  SAiPLE  GPAl’lllC  OUTPUT 


SdOM— MOU 


BEST  AVAIUBIE  COPY 


O 

(O 


o 

CO 


o 

CO 

I 


o 

CD 

I 


CVJ 


CD 


OO 


CO 

QC 


I 

I 

LU 


CD  I 
UJ 


lO 


00 


CsJ 


cn 

CD 

Cl- 

uj 

CO 

CM 


CO 

o 


LU 


O 

cn 

CD 

GO 


O 

(D 


O 


d 

CD 


cn 


o 

CO 

CJ 

UJ 

cr 

CO 

z: 

\ 

f- 

u. 

h- 

>- 

lL 

u.o^ 

UJ 

wOCO 

Q 

-JU- 

UJUJCJ 

a>^ 

3*-*  w 

CD  zoi-i— 

CL— I IU3 

LU  ozctz-i 

ft-)  Ot-iOQt-'U- 

LU 
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INLET  Program  — Subroutine  Control  Structure 


100 


> 


PROGRAM  INLET  ROUTINE  TO  CALCULATE 
EQUIVALENT  K + R FOR  NONPRISMATIC  INLET  CHANNELS 


101 


BESI  AVAIUBIE  COPY 

T *S*  (TNUQ7?a  • TBTS)SI“S 

f TN(R*?) 

ATTiCNfPLO^.PLOTTT, I0e0Q72a«pt“R»l .C7«l ) 

L0*0(PL0T) 

BE«I'JO(TiPE•^) 

COPvSBEfTAPp?, OUTPUT) 

PHOCRa-  I^LtT ( InPuT.OutPUT t TAPE 5»lNPUT.TAPt6«0ijTPuTiTAPt9|T*PEl Of 

1 TAPfc^.PU^CMATAPtT) 

C *<*•»•««•«*««»***«»***«»*«** 

C PO'iTl.f  TO  CALCUIATF  tOUlVALF^T  K ♦ R POP 

c so‘.-ppiS''ATic  Inlet  channels, 

c 

tNTfr.FB  p.  LH,  Pi,  CHlt  P”1 
peal  s.  <•  'P 

CC^-On  inu“''T(lOOO)«0tTAl«niZtTAtAnteAV|OJ»TL*lg»RO*M| 

1 "Z.nT .rTS.Oi AX .mOuBS»E  f5«5)«HlHtOlNOiIOUTiIPLTS»XT 
OI-EnsTON  aC0iH)i  R(A.R).  OLCPtB)*  OOl(9) 

OI'-ENSTfiN  TlTuHu).  tit|,?(u).  EMB.e) 
naTa  JSP/S/ 
haTa  TPLTPP/  i/ 

c 

C the  FCLLn'*ING  ^ CAROS  aPPlT  OnLT  To  CDCfcfcOO  » COPE  1200  CARD  OUTPUT 
C AT  Bf"OTF  Batch  terhj^al 

Punch  1 

j poo“AT r ••»c» ) 
c 

C RPaO  CaBh  type  Ai  containing 

C IP  (I?).  1 IN  CC  ? InOIcaTES  ip  plots  are  OFSIHEO 
C PLANX  INDICATES  No  plots  OESIREO 

read  (ISOlPOOn  IP 
c 

c CH  ■ NO,  OE  Channels 

C P ■ Nn,  OE  X«SECTIONS 

c 

tOUT*», 

IPLTS  ■ 0 

c 

IP  (IP)  10,  10*  s 
s OE"INn  IPLTPE 

CALL  PL0TS( IOUHMT, 1 000. IPLTPE ) 

10  AO  » I,E*2*) 

hRITE  (I0UT,127) 

C 

C READ  CARn  TYPE  A?  CONTAINING 

c CH  (IS).  nu-BER  OE  fLO*  Channels  in  assumed  plow  grid 

c p (IS),  nu-ber  OE  X. sections  selected  across  Inlet, 

C T (»ip,7).  tidal  period  In  hOuRS 

C H (FR_.  tihaL  SE“1. range  In  peet 

c N (EN,(t).  uannIncis  n for  the  Inlet  or  o 
C IF  distinct  values  oP  n EoR  Each  channel 

C And  x-sECTIOn  are  to  he  SPFCIEIEO 

C HAT  (E»,S)«  hTl  surface  area  of  TMf  BAy  IN  SUUARE  FEET  (ENTER  IN 
C E notation  hITh  “anIjSSa  TD  LEPT  AN©  POmER  Of  10  TO  The  RIGHT) 

C peTa  (rH,U)-  DIMENSIONLESS  SURFACE  ARfA  VARIATION  PARAMETER 

C ZETA  (fn,2)-  beach  slope 
C 
C 

REAP  (isnt?002)  CH,  P.  T»  H,  N,  PAY,  RETAt  ZETa 
C 
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t 


ooooooooo 


BEST  AVAIUBLE  COPY 

IPttOFfS)  )<»<>•  J5 

C 0 ^00  ^0,  Of  CMNLS  OB  EOf  b^BK  INOICATFS  fcNO  OF  DAT*. 

C 

15  IFfCH.FO.O)  GO  TO  BB 
C «t*0  C*Bn  TrPf  *3  (?  CAHDS)  COKTaInInC 

C TITli,  TITl2  . (3*i0#A?/3*i0t*2)  - 2 LINES  OF  CMARACTEBS  CAC« 
c output  iofntification,  also  usfo  as  titlf  on  plots, 
c ^OTE•*•  CDC6600  kobo  Sl2F  IS  10  CMARACTfBS 

bead  (150,2003)  TlTLt*TlTL2 
C 

C FOB  P a Oi  OELETF  CALCULATION  FOR  FOUIVaLE'^T  K 

C 

IF  (P.E0,0)  GO  TO  BH 
-BITE  (I0UT,I25)  TITLl,  TITL2 
C BEAD  CaR^  TTPF  Bl  containing 
C LABEL  TO  inpiICaTE  Data  fdLLOkInG  are 
C CROSS-SFCTIONAL  AREAS  (NO  OaTA.  REaO) 

C FOBRATfeOX) 

read  fI50,200A) 

PI  s P ♦ 1 
pxt  a P • 1 
CNt  a CM  ♦ t 
C 

C «FaG/(  1 .BA.ftaaaj 
C 

rF  a lU,5e? 

bead  * PRINT  INPUT  data, 

ARFaS 

REAO  CaRo  type  52  containing 

A(J«I)  • ieF8,0)  • CROSSkSECTIONAL  ARFa  EACM  assumed  FLOR  channel 
(I),  ONE  CaRO  fob  Each  CBOSSkSECT I 0N( J ) , 


00  UO  J a 1,  P 

REAO  (T5««101)  (A(J,n,  I a 1,  CH) 

■BITE  (10UT,lOl)  tA(J,I),  I a 1,  CH) 

AmIn  a 0, 

00  J9  I a 1,  CH 
30  AMlN  ■ A-lN  4 A(J, I) 

IF  (A«tN  • AO)  AN,  un,  <IU 

35  “IN  a .1 

AO  a a“1N 
C 

C retain  smallest  X«SFCTI0N  FOUND  (TmROAT), 

c 

«0  continue 

C BEAD  CaRo  TVPf  83  containing 

C LABEL  FOR  »IDTmS  (no  data  BEAO)  FORmaT(ROX) 

REAP  (T50,200«) 

•RITE  (TOUT, 111) 

■ 0 a p , 

00  aS  J 1 H P 
C 

C aIOTms 

C bfao  Card  type  Bu  containing 

c B(Jti)  ■ (eFs,o)  > ■idths  of  Each  assumed  floh  channel  (!)•  one 
C CARO  FOR  Each  cross«section( J) , 

c 

bead  (ISOtlOl)  (R(J,I),  I a I,  CH) 


oooooooo 


c 

c 

c 

“Z 

c 

c 

c 

uu 

«5 


bo 

C 

C 


WRITE  flOUTilOl)  (R(JiI)»  I ■ (•  CM) 

00  ttb  1 ■ I#  CM 

TF  (J  . “IM)  itU,  itf,  itu 

COMRUTF  msl  wIOTM  at  ImbOaT 

wO  ■ mO  a R ( J » 1 ) 

convert  «inTMs  TO  Mvnw,  Raoii, 


best.  AVAILABLE  COPV 


R(J,I)  ■ AfJiI)/R(J,l) 

CONT  I»J|IE 

00  so  J ■ 1 I R1 t R 

00  50  I ■ 1 f CM 
OL(Jtl)  ■ 0, 

CONT  TmUE 
TL  ■ 0. 

REAO  CARn  TYPE  Bb  CONTAINING 

label  for  lengths  (NO  Oat*  re*o)  FoRMArreox) 

Rf*0  (T50.200<t) 

00  52  J ■ 1i  PM) 


lengths 

REAO  CaRo  type  B6  containing 

OOLd)  (8F8,0)-  lengths  between  C'^OSSwSFCTIOnS  along  the 

HOliNOARIES  of  The  ASSUMfO  flow  CH*NNfLS(n,  the  number  of 
lengths  per  C*R0  is  one  greater  than  the  number  of  channels 
ANp  the  number  of  Cards  is  one  less  than  the  number  of  sections 
READ  (150*101)  (OOLd),  1 ■ 1.  CHl) 
oo  5?  I » 1*  CH 

OL(v'*t,n  1 0.b»(D0L(Il  ♦ OOL(lFl)) 
yl  » Yi  ♦ OL(J»i,n 
52  T*  ■ I 

n^uYfUS) 

I'U  5'i  J • 1 • Ml 

5i4  hNITE  (lOuT.lOl)  (OL(J,I)»  1 ■ 1*  CH) 

TL  ■ Tl/CH 

c 

c FOR  N 7FR0*  READ  DIFFERENT  N'S  FOR  E*CH  *»SECTI0N, 

C 

IF  (N)  l5b«  156*  ISa 
lb«  90  155  1 P 1*  CH 
00  155  J « I*  P 
EN(J*I)  ■ N 

155  continue 
GO  TO  15* 

C REAO  CaRo  type  B7  containing 

C label  for  maMnGS  n (NO  data  READ)  FORMAT(eOX) 

c NOT  used  if  value  of  n on  c*Ro  Type  *2  is  not  zero, 

156  READ  (I50.200R) 

•RlTF  (I0UT*121) 

00  157  J » 1*  P 

C REao  CAPn  TYPE  H8  containing 

c EN  (J.ll  • mannings  N For  each  channel  of  each  CROSSwSECTION, 

C NOT  usto  IF  VALUE  OF  N ON  CARO  TvPt  *2  IS  NOT  ZERO 
REAO  (150,121)  (EN(J*I1,  I ■ 1*  CM) 

xRlTE  rinuT,l21)  (ENIJ,)).  1 ■ 1*  CH) 

15T  continue 

15*  ROOT  > 0, 

00  60  I ■ I . CH 
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BEST  AVAILABLE  COPY 


ROOTI  a 0, 

00  bS  J » 1 I P 

ROOTI  a ROOTI  ♦ tDL(J*I)  ♦ DL ( J* 1 , IJ ) aR ( J , I ) * * ( • 1 , 3 J J J? 3 ) 

1*(Em( Jf I)/A( Ji  I)  )**? 

5S  COMTiNUfc 

ROOTI  a l./(S0PT(J,/f*(MlN,n»a2)  ♦ KraHOOTI)) 

ROOT  , root  ♦ ROOTI 
60  CONTIkiUt 

c 

C calculate  equivalent  k 

c 

K S 57?i'>6*T*H00T*S0RT(6«.3*M)/(H*HAr) 

C 

C CAI.CUlaTE  equivalent  R 

c 

RO  a 1 JOS.RTTaAO/IKAHAV) 

RO  s ih,O7S*P0*t<n/M  > 1, 

RO  a fn,niu33»*R0/(N*N*TLn»*(»0,75) 

C 

C OUTPUT  data  ♦ RESULTS 

C 

-RITE  (IPUT,103)  T,  H,  bay,  N 
"RITE  finuT.loS)  K,  AO,  Tl,  ROi  RO 
C CALL  ROuTINiE  to  solve  EQUATION  SET 

CALL  I“LT2tTITLl,TITLdI 
GO  TO  10 

C REAO  CaRo  type  Au  containing 

C AO  ffS.O)-  mTL  CwnSS"SE CTIOnAL  inlet  area  in  square  feet 

C .O  (TS.O)-  TNLET  W-IOTM  I».  EFET 

C RO  (TR.a)-  INLE^  OEPTm  or  hyoWaULIC  RAOTUS  IN  EEET 

C TL  (AS.O)-  INLt^  LENQTm  in  feet 

R«  read  ft^jOiPOOS)  AO,  «U,  «0,  Tl 

CALL  I*'LT^fTITLt.TITL^l 
GO  TP  10 

R9  IK IPlTS,CT,0)  call  PL0T(XT,0,,R99) 

C 

C TRE  EOLLP-ING  2 CARDS  aPPlY  ONLY  TO  CPC66flO  • COPp  1200  CARO  OUTPUT 
C AT  remote  Batch  TfRI'Inal 

PUNCH  QR9 

999  FO«HAT(a«PL*) 

STOP 

c 

101  FORhAt  (BEP,0) 

103  E0R“At  (/  ISH  TIOal  PtRlOD  a,  E7,3,  3rt  mr/ 

llBH  TjOAL  SLhiRangE  a,  E7,3,  3H  Et/ 

217H  B«v  AREA  (HSl)  a,  1PE12,6,  6H  SO  FT/ 

3mH  MANNING'S  N 3,  0PE7,«/) 

C 

105  FORMAT  nSH  COEE  OF  REPL  «,  F7,3/ 

lliH  Inlet  area  a,  1PE1?,6,  6H  SO  FT/ 

215h  Inlet  length  a,  Fl?.h,  3h  ft/ 

SIRH  I,, let  oEPTH  a,  0PE7,3,  3h  FT/ 

<I1UH  I»lET  mIDTh  a,  F7,o,  3h  et/J 
C 

lOT  format  (5(lPEll,S,  1H,1/  ifE11,b,  1H,),  FH.S////) 

111  format  (/7h  mIDTmsJ 

US  format  (/ew  LENGTfiS) 

121  format  (/  lUH  MANNING'S  NIS) 

123  format  (BFB.D) 

125  format  (/5X,  35mEQU1VA(.ENT  COEFFICTENt  OF  REPLETION/ 

lll*t  23HFaR  NON-PhTSmATIC  IMET  // 
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BEST  AVAIWBIE  COPY 


2 2f6X,  3*lf'l*2/)/  hK  ABfcAS) 

c 

J27  FOR»-at  (/UOMI*  •••*•«*»«•«»««*»***  /) 

2001  FGH^at  (12) 

2002  FCRmat  {2I*it  MO. 7.  Fb.Jt  Ft>,<i.  F8.3«  F6.«t  Fo,2) 

2003  FOSf‘AT(3Ain.A2/3A10.A2) 

200u  KiUmat  (MOX) 

2005  FOS^-AT  (2F8,0.  M,Ki  Ffl.O) 

ENO 

SllBRniiTINf  InuT2(ai,A2) 

c 

c «**•« 

C P»0r,B4“  I^LtT^  » AHAL^StS  T^lt7-flAV  SrsTfK 

C RV  l."i“PEr)  PaHA'^LTEP  APPWOaCM.  I*»CIU0E3 
C InEPTja  of  FLOA.  ALL(I''S  FRfSh  «aTfb  InfuDa, 

C 

c oeeinf  ocEA^  TIDE  AND  tnfloa  as  Functions  of  ti^e  in 

C SUBHOilTINt  Ff»T(T.M2) 

C T ■ TiA't  IN  SEC, 

c Mi  ■ ocean  tide  in  feet 

C 

peal  xF.  n,  L 

dimension  41(u),  a2(u),  DUMMV(b) 

COMMD>  IOi.'““T(10nn),RET4*A0,ZETAt4X0*AR0.GLtLfN»P0« 

)M,XF.oT,OT5,fJ,A«,HDLOS,F(5,5).MlM,OINo,IFLlfIPLTS. 

2»T 

DaTa  lFLi/S,7/,  1 ,F  ♦ 75/ 

C 

C HEAD  CaHd  type  A5  containing 

c T (Fr,5).  time  at  amicm  calculation  IS  TNIIIaLIZEO 

C Tn  Fh.si-  time  at  ahIcm  piJTPUT  IS  TO  HFGIN 

C TF  Fr,5).  time  In  hOI'HS  aT  »mIcm  SOLUTion  IS  TO  BE . T E Hm i n A T FO 

C Mjo  (Fb.3)-  initial  condition  fiiH  BAY  surface  ELEVATION 

0 V (Fb.3)-  initial  condition  for  inlet  velocity 
C DT  (Fb.o)-  IImE  step  In  mInuTFS  To  he  used  IN  NUMERICAL 

C SOLUTION  method 

C NOT  (In)-  specifies  Nhmh(.R  of  Tlnfs  STfPS  TO  BE  PRINTED 
C IPLT  (!«)-  SET  EUUAL  to  1 IF  SOLUTION  PLOT  IS  DESIRED 
read  (IFL2.2010)  T* TOfTF •Ml 0» VtDT.NDT , TPLT 
IF  (IpLT.nE.O)  HE*IND  9 
mRITE  (IFl1»100U)  ai.  a2 
rHITE  (IELt.1005)  M,  arOi  beta 
IZETA  c ZETA 

RRl'TF^dFLt.lOOb)  AXOf  aO,  IZETA.  PO*  Lt  N 
aHITE  (IEL1*1009) 

DT  = dO.AOT 
GL  B J?,15/L 
RE  a ?9.122«n*N»l 
beta  a RETA/m 
T s 3600. »T 
TO  t 3600, *T0 
TF  a 3600, *Tf 
DT5  a o.SBDT 

CALL  DEO2(T.Hl0,Ml 1 ,v*vl IHJO) 

ROUnT  a 0 

230  IF  (T  - TO  ♦ DT  ♦ DT5)  2“0,  ?uOi  2S0 
2«0  CALL  OEO?HR(TfH10tHll*VfVl,M20) 

ROUNT  a ROUNT  * I 

C 

C safety  valve  • ABORT  IF  ROUNT  EXCEEDS  50OO 
C 
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BEST  AVAIUBLE  COPY 


IF  (Kni)M  - Sooo)  iiOt  2J0.  2uS 
2«5  WHITE  (IFli.1007) 

HETuHnj 

C 

C BEGIN  SOLUTION  CYCLE  "ITM  OUTPUT 

C 

250  E{  1 , 1 ) 3 H?0 
E( I t?)  r OINO 
EIl.J)  * M10 

EClfU)  3 V 
FdtSi  3 0 

CALL  OFu2HK(T|Hl0,Hn«V|Vl,M?0) 

Eta.lJ  3 H?0 

F(2.2j  3 tJiNO 
F(2|J)  » HlO 
E(2iU)  3 V 
Ft?, 5)  s 0 
W1SU“  3 0, 

CALL  OEU2K'<(T,H10,wn  »V,Vj,m2o) 
louT  3 n 

CALL  CHIT  tM20,Hl  o.V,  IPWNT,  IOL'T) 

IF  (IPLT)  ?71,  271,  256 

256  WHITE  flELI)  hOuHS,  tEn,lA),  lA  s 1,  5) 

GO  TO  271 

260  no  27n  I ■ l»  NOT 
lOIIT  3 0 

Hisu«  3 wisuM  ♦ Hin 

CALL  nEG2HKt T ,M1 n,Hi 1 ,v, VI ,h2o5 

CALL  CHIT  f h20,Hl 0, V, IPHNT, lOUT) 

IF  (IpLT)  270,  270,  265 

265  WHITE  (IELT)  mOuHS,  fFfl,U),  U m ],  5) 

270  CONTInIpE 

IF  (IniPT)  271  , 271  , 260 

271  WHITE  (IFL1,1010)  houhs,  (E(1,IA),  IA  s i,  5J 
280  IF  (T  - TF  - OT)  260,  2<>0,  2H0 

290  HlSU”  3 hi8UN»0T/(TF  - TO) 

IF  tiPLT.GT.O)  white  IIFLJI  OlJMMY 
wHITF  fIFL1,10in  wiSUM 
C 

C END  Of  solution  routine,  on  TO  PLOT  HOUTINt  OH 
C BEAD  next  data  set, 

C 

IF  CIPLT.GT.O)  call  GHPmC(31,*2) 

RETUHo 

c 

1001  FOHmat  (3(5A6/)) 

lOOil  format  (IHl*  20X, IRmINLE I ANAlYSIS//?(12X,«a6/)/) 

1005  format  (5X,  17HTI0AL  SfMlhANGE  3,  F7,J,  iH  FT/ 

15X,  23WMEAN  BAY  SURFACE  AREA  s,  IPFIO.R,  6H  SO  FT/ 
25*»  31HBAY  SIDE  SLOPE  PaHam£TeH  beta  3,  OPF6,i/) 

C 

1006  format  (5X,  17MINLFT  PhOPEHTIES|/  BX, 
l?6F'X-SFCTinN  area  HELCw  mt|  3,  lPEll.5,  6P  SO  fT/ 
2«X,2nMlNLET  wIOTm  at  m.tl  3,  PPFT.o,  3m  FT/  8X, 
ilHBINLET  BEACH  SLOPE  «,  lu,  2miI/  ba, 

UIIHMTl  depth  3,  F7,3,  3H  FT/Bx,  iRHlNuET  LEnGIh  b, 
5F7,0,  3H  FT/  dX,  ISHMANmIngiS  n b,  FT.u//) 

C 

1007  format  (29m  execution  aBOHTEP,  TRANSIENT/ 

126M  stage  exceeds  5000  STEPS,/) 

C 
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BESOVAllABLE  COPY 


1009  F09m*t  («*•  UhTIMti  S*.  ?MHa.  UXt  fenKFLOKt  bx,  aHMl, 
1SX«  3hVEL>  7X«  IMQ/  Sx,  3nHWS,  bXt  ZHFTi  5Xi  UmKCFSi 
16X|  2HFTi  sx.  3hFPS,  6x,  (imkCFS/) 

C 

1010  Foat-XT  (2FF,at  »3aF9.c!,  0P?F«,2,  -3PFin,?,  ^^9,3) 

1011  F0af**T  (/  lb"  *vc  RAV  LtVtL  «.  Fh,?.  3"  FT// 

123"  * CRITICAL  Pni>JT  VALUE//  lOXt 

228"*A*  TABULATION'  CU"PLETtn  *••/) 

C 

2010  format  (3FR.5*  2Fb.3.  Fb.O.  Ib«  Ul 
ENO 

SURROuTINE  0ES2  ( Tt"l 0,M1 I , V, VI ."21 

C 

C SURRnllTl^E  UtU2  ASSEMbltS  AND  EVALUATES  T"E  N0n»LIREAR 

C simultaneous  OiE.'S  C"ARACTERi/Inc  T"E  INutT 

C 

real  kF,  Lt  N 

common  inu“’'T(l000)«RETA»«OiZETA»AX0iAR0#GL*LtNfR0*"iKF» 

C 

10T«OTs,OtAX.HOURSiF(3»S).MlM,olNOtTFLl.IPLTS*XT 
C CALL  FEXT  FOR  tide  AND  InFiO* 

C 

CALL  FfXT(T,M2l 
ETAO  e O.S«("lO  ♦ "25 
C 

C VELOCITY  COEFFICIENT 

C 

CV  » n,01bS5*(l,  ♦ xF*(RO»ct«o)*»(.1 .33333) ) 

C 

c Surface  area  of  bay 

c 

A8  « aR0» ( 1 , ♦ be TA»"10) 

c 

C INLET  cross-sectional  AREa 

C 

AX  ■ AXO  ♦ (-0  ♦ ZF TAAfTA0)*ETA0 
VI  • r.L»("?  • "10  • cv*v»ARS(v)) 

C COMPUTE  OV/OT  Fhom  ENERGY  FOUaTIOn 

C 

0 • A**V 

c 

C compute  Omi/OT  FRO"  CONTINUITY  EQUATION 

C 

“11  ■ f Q ♦ GINO ) /AB 

"UoRS  ■ (2,77  )7777E-0IiI*(T  • DT) 

"1"  ■ 3,bE*0b»"ll 

PE  Turn 

FNO 

subroutine  FEXT(T.m2) 

C 

C subroutine  to  define  ocean  tide  "2  aNq  INFLOm  OINO  AS 
C functions  of  time  t in  seconds 
c 

real  kF,  l.  n 

common  I0UM"V( 1000.) iBETa.-O.ZETa. AXO. ARO.GL.L.N.RO.m.KF. 
IOT.OT5.O. AX.MOURS.F (3. S) ."IM.OINO, IFLl .IPLTSiXT 

c 

C CONVERT  SECONDS  TO  RADIANS 

C 

X ■ 1 .<I0S63RE>0««T 

c 
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C SINIjSnID*L  not  OF  *HPUITuPt  H 

c 

H2  ■ h*sIn(X) 

C 

C IKFlO.  sEGLiGIBLt 

C 

OINO  « 0 

BE  TuS>' 

suBunuTiNf  cpit(h20»kip,v,tbbnt, loiTr ) 

c 

C SURPO'iTIMt  cwn  CO''P*'<FS  3 CflMSECUTIVf  FMNCTIOm  POIi-TS 

C »NO  -BTTtS  •'IDDUE  POI^T  IF  II  IS  A CPTTICAL  POINT 

C 

real  kF,  Lt  N 

Cll^MOtj  IOU»'NV(lonO)fB£TAt"PtEETA«AXO*AnOiGL.LtN#RO|HtFF| 

IPT.OTs.O. AX.HOURStF { 3.5) iHlH.OlNO, IFUl , IPLTSf XT 
niNtNSION  «aRx(5) 

OATa  uARXA/IH  /,  MARKb/lH»/ 

F ( 3, n a M?n 
F(3.?l  8 QInO 
F(3.3)  8 HJO 
F( 5.U)  a V 
F(3.5)  3 0 
00  ?n?0  TA  8 1.  s 
HAR»(TA)  8 MAWKA 

IF  (F(?,IA)  • F(t,lA))  201?,  2020.  20la 

20l2  IF  fF(3.lA)  • F(2.IA))  2020.  2Pl5,  2015 

20l«  IF  (F(3,Ia)  • F(?.Ia))  2015.  2015.  2020 

2015  TOUT  a 1 

«AHk(ia)  a MANKH 
2020  CO*ni>..tit 

00  202F  lA  a 1.  5 
FII.Ia)  s Ft2.lA) 

2025  F(?.Ia)  a F13.1A) 

TF  ( Tout  ,E0, 0 ) HFTliRN 

.RITE  (IFlI.2101)  hoiiRS  , ( F ( 1 , I A ) . m aBK  { j * i , I Aa  1 , 5 ) 

BtTUH*J 

2l0l  FciRPAt  (2F».2,  Ai  ,.JpF8.2.  A|  .2(0pF7,2.A1  ) . 

4-3PFB.2.  Ai,  2(FP,3.  *1)) 

END 

SUBROtiTiNE  DFU2RK(X.T.V1.2.Z1.S) 

C 

C RUNGE-fuTTa-CILL  SiirbOiiTInE  to  solve  2N0  ORDER  O.E. 

C OR  S,  sTpulTaNEOUS  first  OROtP  O.E.'S 

c 

real  aF,  L.  N 

common  lOlJMHydOOOI.HETA.afl.ZFTA.AXO.ABO.CU.U.N.RO.H.KFl 
10T.nT5.Q. AX.MQUHS.F' 3.5) .HIH.OINO. IFLl .IPUTS.XT 
CALL  riFO?(x.T»Yi.Z.Zi.S) 

PI  a v1*0T5 

01  8 218015 

X a X 8 0T5 
YOl  a V ♦ PI 

Z01  8 7 ♦ 01 

CALL  PFt)2CX.Y01.Yl.Z0J.Zl.S) 

P2  a V18DT5 

02  a Z180T5 

YOl  a Y ♦ P2 

ZOl  a Z « 02 

CALL  nEQ2(X.Y0l .Yi .ZOl .Z1 .S) 


P3  • yt*OTS 
Q3  « Z1*0T*i 
X ■ X ♦ OTS 
vci  ■ V ♦ P3 
ZOl  • Z ♦ f33 
CALL  nF.U?(X,YOI|Vl.Z01,ZliS) 

Pu  ■ v^priT? 

QU  B 7i*l'T5 

ni  B f>.3333i333*fPlAP«i*P2»P2*P3AP3) 
f'2  ■ n.i3i33333*(0U04i*Q?*O2A03*O31 

V B V ♦ ni 

Z » Z ♦ 0? 

CALI  nK)?(x.v,yi ,z.zi *8) 

PtTuSn 
* NO 


SUPROllTISl  PLOYLNtX.V.l.N*  IMOX) 

c 

C SiiRHOiiTiMt  TO  Plot  vapiohs  lImE  types 

c 

IPEn  a 2 

TNOX  B INOX  ♦ 1 

GO  TO  (20tP0«30<10)f  LN 

C .•••••«  L^iPP 

10  GO  TO  (uOt  1S)|  INOX 

IS  T*yOX  a 0 

GO  TO  3S 

C ...  . ...  . ...  LN  • I 

2o  GO  TO  r«n,u0tU0t3SiP0t2S)«  I^OX 

25  iNOX  B 0 

GO  TO  35 

C ...  ...  ...  ...  LN  • 3 

3o  GO  TO  run«35tU0t35tP0«35iPOf P0»«0«32) t INOX 

32  INOX  B 0 

35  IPtN  B 3 

C ...............  LN  « 2 

PO  CALL  PLOT  {XtYlIPtN) 

return 


END 

SUHROuTINE  AXlSt3(XtYi IHCDiNCfHiSIZEiNNt  TXVtXMlNtOXiSPACEt iTICt 
INPIO) 

niNENsION  LHL(3) •NfSAGPCSj 

oaTa  lPL/Th(X10  )/  1 

OATa  “ESAGE/liNAXIS  TOO  LONG/ 

XYSIXy 

AnGL  FbXYARO . 

TICsItTC 

SPACbsPaCE 

r.sM 

PIObNpJO 
XTBX 
yTby 
lINf  B? 

SIZbSIZE 
IKST?)  I f IN*  loo 
1 LInEbt 
SIZb.sIZE 
loo  nQbnA-nn 

IKNO,  J,2,? 

2 NDlGBNP** 
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BK  AVMIE  COPY 


r.o  TO  ^0 
j ^^n•^nTr.■o 
50  sS«>*Cf»SW/SP*C*,5 
Tl.»>.SP»Cf  •SP*C 
tFfYT*Tu»«v-i<».) 

50  C*LL  S»‘'5OL(*TiYT*,i,,i«,mfs*CE.90,i15) 

QtTo<»>, 

<1  •<C?sr.»,5 

H07T?u««i»r. 

PO^t  9s  1 0 t»»  t •'*910  ) 
r5tL»iP'7»E«»OX 

f(apn«FS*XMIN 

*‘.G«  1 .-XY 

MbPNC 
*L*H*l  . 

17  ('-BjSiO.O 

5 NBSaNR 

6 TIC»rsf2,»*Y-l.)**L*B*TlC*,l 
*TIC*TTCXV*XY 

YTIC»T1CXV»***& 

ANSNiH 

17 (Pln)60.65»60 

bO  XKSXK^O, 

65  8TnLFeTL».5-(7,»4K.3,l»H07*,5 
TIC'"1a.lh-TIC*.05 
4*iO“I>,«2U, 

XM)'-Rb)(''E  •* 

NMXXBfl 

I9X«fl 

00  I?  I*Y9a.mS9xCE 
NOlGlTB^Dir. 

CALL  PI 0T{XT.YT,LIN) 

Ll'*aLI‘*E 

CALL  Pi OTCxTYXTIC.YTaVTICi?) 

17(G)  7t 1 1 ,7 

7 ITE'-Oi,<45u?9U«e2*AL06lAHS(XMJMH)*,5*J0,**{«ND)  )♦!  , 

17 (ITfvp)75t75»5l 

51  COMIsUE 
noigI7«'ydigit*ite"9 

75  T7(x'*i.>'H)9,8tlO 
e TiDiGiTtO 
9 NDlGlT*NOir.lTYl 

10  I7(f.0lGIT-LMAX)  77*77.5? 

52  CO'*Ttk)llt 

|gMAxa^r)IGIT 

77  CEMTEPsM07»(7,*NniClT-3,)*,5 

XAMO*-AL*H*xY* (T ICKlACfNTEPl-cf Mer 

VA'*0«aLABAANG*  (TICWlAbnpl-MOa 

CALL  number (»TYXANO.¥Y*YASO,G,*NUMR.O,,N4) 

XNUMBxXNUmraOELX 

tF(ANO“l'*-*ANO)  11*11*55 
55  anO»In»xano 

11  CALL  PL0T(xT*yT*5) 

XT«XTaSPAC»AN6 

YT«YT*SPAC*XY 

12  COMTIUUE 
T7(G)i5*16*15 

13  ANOi»Id«(7,9NMAX»3,)»mO7*0,5 


III 
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BESTAVAIUBIE  COPY 

)(TITLt*XVYV»(*N0MlN+AN0t«lD-*l.*H*(*N0WID*,l8))»*NC*8TITLt 
XT ITLE»X**MG* (XAM0+ALA8* (CA, 1 6) )+X**ST tTLt 
CAUL  syMhOL(XTITLt«VTlTLt,GflBCn*ANGLtfMB) 

IF{PIOU««16»l« 

\U  BM(Nb»1)»G 

XT  jTI,f«XTITLEtANG*HN 
XTlTLESYTlTUtAXYAB^ 

CALL  SVMBOLCXTITLEtTTITLt .G.LHLC 1 )# ANGLE. 7) 
XTlTLF«XTITLEAANG*H07T2u-XV*Mn2 
vTITLt«VTITLE  ♦XV»M07T?aAANG*MOa 
CALL  NMM8tR(XTITLE.TTlTLE.B,*H07.Pl O.ANGLE.-I ) 

Ifc  PETuPn 
END 

SUBMOiiTINE  kEAOIn  ( x . Y , y E AC  . X E AC  . X 0 • XE  , I mOC  . KK , LN) 

C 

C subpOmtine  to  reao  solution  Tabulation  from  file 
C 

real  kF,  l.  n 

COnmDn  inuMNY ( 1 000) . beta.  »0.  ZETA.Axo.aBO.CL.L.N.HO.H.KE. 
lOT.OTS.Q. AX.houRS.E ( T.S) .MIH.OINO. IFLI . IPLTS.XT 
OINENSION  Y(5).  YFAC(b) 
read  (Q)  X.  Y 
IE  fX.GE.I.Et/S)  KK  ■ ? 

INOC  a 0 

IE  (Kk  - 1)  lO.  10.  SO 
lO  IE  (Xn  • X - OTi)  20.  SO.  SO 

20  IF  fx  - XF  • OTS)  10.  25.  25 

25  KK  a ? 

GO  TO  SO 
JO  INOC  a I 

X a XFAC*(X  • XO) 

Y(LN)  a YFAC{LN)*Y(LN) 

So  return 

END 

SUBROUTINE  CRPMCCALABLI .ALAB12) 

C 

c subroutine  grphc  -rites  Plotter  tape  for  GRAPhicAu 

C output  of  solution 

c 

CO“nOn  I(VU“MY(  1 000)  .RUliMBYf  12)  .dT5.S0ummy(20)  . IEL1.IPLTS.XT 
OI“ENSION  YLAbLLEJ) 

DINFNSION  alEGN('i.6).  aLABI  ( u)  . aLARL2(<i) 

OI«ENSTON  SYN(l),  Y(S)..  YFaC(5) 

oat*  xLABL/  RMTlMt-«HRs/«Y|  AHLR/10nEL0'.--XCFS/ 

OAT*  yLABLU/IOMneIGhTS,  V , l owe LOC I T I E S« . BH-F T . EPS/ 

OAT*  *LEGN/10hFLO»  (kCFS.IOM)  ,Jm  .IOMINLET  VELO.IOhcITY 

1 (ET/S.IMEC) • IOmBaY  TIoE  (.IOmET)  .1M  .IOmInFLO"  .lOH 

2 .IN  .IOMOCEAn  TlOE.lOH  (FT)  ,iM  .lOMLECENO  .lOH 

J .JH  / 

TPLTS  a IPUTS  ♦ I 
OTS  a OT5/Jb00. 

C HEAD  Caro  typf  ao  cont*ining 

c XO  (FB,5).  starting  TIMF  OF  PLOT 

C XF  (FB.S)-  FnDING  TiHf  OF  PLOT 

C SCALX  (FB.S)-  TI^E-AXIS  SCalE  IN  mOuRS  PER  INCH 

c YLO  (6,1)-  PINIHUR  VALUE  (lE  tIOaL  HEIGHTS  IN  feet  OR  VELOCITIES 

c IN  feet  per  SECONO 

c YL  (F6,1)  - overall  height  of  plot  TN  inches 

c YLSCal  (F6,1)-  scale  of  tidal  heights  in  feet  PfH  inch  or 

C velocities  in  feet  per  second 

c VRO  (Ffc.ll-  HINIHUH  VALUE  OF  FlOmS  IN  THOUSAND  CUBIC 
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BEST  AVAILABLE  COPY 


C FEtT  PEW  SFCONn 

C YHSCaL  • SCALE  OF  FLOW  IN  TmouSawP  CUBIC  FEET  Pt»  SECOND 

C SCALE  (P6,J)  - scale  EACTOW  EOr  TOTAL  PLOT  SIZE 
PEAD  ( S<20Ol)  XOf XFtSCALXf TLOtYLf VLSCALtYROi 
JYHSCALf  scale 
CALL  r acTOP(SCALE) 

CALL  PLOT(O,i0tt»i) 

C 

C PLOT  legend 

c 

no  2P  LN  » J*  j 
iNnx  B 0 

YP  ■ O.B  ♦ 0,2»LN 
CALL  PLOT  f0.»YP*0.06*S) 
no  IS  I ■ It  p 
xw  ■ 0.1*1 

I INt  ■ « - LN 

15  CALL  PI.OTlnCxP.  yp+o,ON.  LINtt  INDX) 

SYM(l)  > ALtGMl.LN) 

SYW(?)  B ALEGN(2.LN) 

SYM(Jib  ALEGN(J.Ln) 

CALL  SYMBOL  (l,liYP.0,l(ItSYM,0,,2J) 

20  continue 

C 

C print,  1ST  20' 

c 

MHlII  fIFLltl0l)2) 

no  in  LN  ■ «*  *> 

|NDX  • 0 

YP  e O.B  ♦ 0,2*LN 
CALL  PLOT  f0..YP*0.0Sii) 

00  25  I ■ 1 I P 
XP  ■ 0.1*1 

25  CALL  PLOTLN(XP.YPF0.OB,aiINDX) 

SYM(l)  s ALEGNd.LN) 

SYM(2)  a ALEGN(2tLN) 

SYW(i)i  AlEGN(5.LN) 

So  CALL  SYMBOL(l.liYPtO.UiSYM|0.t23) 

C 

C print,  1ST  iOl 

C 

WHITE  flELtilOOi) 

CALL  symbol  (0.6t2.1t0,21tALEGN(ti6)t0,i6} 

CALL  symbol  (3.5f l,75tO,21,ALABLlin,f3?) 

YMIO  B i.O  ♦ 0,5*YL 

TOP  ■ 3,  ♦ YL 

XL  • (XF  - XO)/SCALX 

CALL  symbol (“,25t I ,25i0,2l . ALABL2.0, t3?) 

C 

C PRINT,  'PLOT  X AXIS' 

C 

WRITE  flELlilOOP) 

XTIC  B 1,/SCALX 

XL13  a -XIIC*ELOAT(IFIx{XF  • XO)) 

53  CALL  AXISli  (0.i3,0tXLABLi«9t0.1O|XLl3t>tt0tl(0, 

It.OtXTlCtliO) 

C 

C PRINT,  'PLOT  NIGHT  AXIS' 

C 

WRITE  lIFLltlUOS) 

OYR  a #0.5*790 


BESUWIIABIE  COPY 

> i*  • . • > ' • ’ 


YHTIC  ■ nYP/rHSC*L 

CALL  AYISM(YL>i.O(YLARLHflC|0.|lliYLfliltYRO| 
iriYB,  vBTICiItO) 

CALL  PLOT  (ALiTOP,^) 

CALL  PI UT  fO.*TOP,a) 

C 

C PRIM,  iPLOT  LEET  tXiSl 

c 

OWITF  (irutlOOt,) 

DYL  • -OtSYYLO 
YLTIC  ■ OYL/YLSCAL 

CALL  tXISM  (O.tJ.n.YLABLLtZStO.KliYLfttl'TLO, 
lOYL.  Yl  TTC.  l»  0) 

CALL  PLOT  fO.tl.OiS) 

CALL  PI OT  (XLti.Pt?) 

CALL  PIOT  (ALtYxiDfi) 

CALL  PLOT  (O.fYMin,.?) 

c 

C PHiNT,  iPOBOfW  COMPLETr* 

c 

call  pi 0T(-.?^«0,*»1) 
wWI TF  f IFL 1 • 1 nOT ) 

YFACd)  • l./YLSCAL 
YFACC?!  ■ n.UOl/YHSCAL 
YF  AC  ( SI  ■ YF  AC  ( n 
YFACfu)  ■ YFAC(l) 

YFACfS)  • YFACCiJ 

IMF  « (I 

XF  AC  "■  I ,/SCALX 
C 

C PLOT  nctAM  tide  ANn  imflOi* 

c 

OU  BO  I ■ li  j 

c 

C PWlMT,  I no  BO  1 ■ • t I 

c 

xRTTE  flFLl.lOOR)  I 
CALL  PLOT 

KK  S ) 

iMOx  ■ n 

35  CALL  RFAniM(X«YtVFACiXFxCiXOfXF,lNnC«KKtn 

CO  TO  filO#P5)i 
«0  IF  (Unc)  SBi  35*  P| 

R1  call  PLOTlMX.TCD.k.ImdX) 

IF  {Y(2))  <i?.  35.  ui 
«2  IMF  ■ I 

GO  TO  SB 

<•5  IF  (InF)  60*  60*  50 

So  C0BTI(.;llE 

C 

c Plot  ray  rfspomse*  velocity*  flo" 

c 

60  00  eo  I ■ 1*  3 

c 

C PRIf.T*  100  60  I «'*  I 

c 

»»ITF  tlFLl.lOIOJ  I 
CALL  PLOT  (U.*  U.*  3) 

g 
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BEST  AVAILABLE  COPY 

iNDx  z n 

6S  CALL  RM1IVI  (XtYtYFACiXf  *C.XO,X^*INOC,KK,I»£) 
r.O  Tf)  f70,  'eo).  KK 
70  IF  (I‘inci  hSt  65*  72 
12  CALL  PL('TLM(X|YCt*2).I,lNDX) 

GO  TP  6b 
80  CO'JTtMif 

C 

C PH!N7,  1ST  60' 

C 

• *<m  (IFLIflOOB) 

XT  » XL*". 

call  PLOT(xT.0t»»3) 

PETUB^ 

c 

1001  F0»'‘AT  (6A6) 

1002  format  (6M1ST  20) 

1003  format  (6m  ST  30) 

lOOu  FO.Rmat  (12m  FlOT  X AXIS) 
loob  format  (16M  Plot  rigmt  axis) 
lOOb  format  (ISM  PLOT  LFFT  AXIS) 

1007  format  (IbM  HORDfW  COMPLETE) 

1006  format  (6M  ST  80) 

lOOR  format  (lOM  00  50  I Bf  12) 

1010  format  (lOM  00  80  I 12) 

200l  FORMAT  (3F8,5.  6F6,3) 

EnO 

EOR 

COR 

COI 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Huval,  Carl  John 

Comparison  of  numerical  and  physical  hydraulic  models.  Mason- 
boro  Inlet,  North  Carolina;  Appendix  4:  Simplified  numerical 
(lumped  parameter)  simulation  / by  C.  J.  Huval  and  G.  L.  Winter 
gerst.  Vicksburg,  Miss.  : U.  S.  Waterways  Experiment  Station  ; 
Springfield,  Va.  : Available  from  National  Technical  Informatioi 
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